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ABSTRACT 

We investigate the use of mid-infrared (MIR) polycychc aromatic hydrocarbon (PAH) bands, con- 
tinuum and emission hues as probes of star-formation and active galactic nuclei (AGN) activity in a 
sample of 100 'normal' and local {z ~ 0.1) emission-line galaxies. The MIR spectra were obtained 
with the Spitzer Space Telescope Infrared Spectrograph (IRS) as part of the 5pjteer-SDSS-GALEX 
Spectroscopic Survey (SSGSS) which includes multi-wavelength photometry from the ultraviolet to 
the far-infrared and optical spectroscopy. The continuum and features were extracted using PAHFIT 
(Smith et al. 2007), a decomposition code which we find to yield PAH equivalent widths up to ^ 30 
times larger than the commonly used spline methods. Despite the lack of extreme objects in our 
sample (such as strong AGNs, low metallicity galaxies or ULIRGs), we find significant variations in 
PAH, continuum and emission line properties and systematic trends between these MIR properties 
and optically derived physical properties such as age, metallicity and radiation field hardness. We 
revisit the diagnostic diagram relating PAH equivalent widths and [Ne ii]12.8/xm/[0 iv]25.9/^m line 
ratios and find it to be in much better agreement with the standard optical star- formation/ AGN clas- 
sification than when spline decompositions are used, while also potentially revealing obscured AGNs. 
The luminosity of individual PAH components, of the continuum, and with poorer statistics, of the 
neon emission lines and molecular hydrogen lines, are found to be tightly correlated to the total in- 
frared luminosity, making individual MIR components good gauges of the total dust emission in SF 
galaxies. Like the total infrared luminosity, these individual components can be used to estimate dust 
attenuation in the UV and in Ha lines based on energy balance arguments. We also propose average 
scaling relations between these components and dust corrected. Ha derived star-formation rates. 
Subject headings: surveys ~ infrared: galaxies - galaxies: star formation, active, ISM 



1. INTRODUCTION 

Determining the main source of ionizing radiation and 
the star-formation rate (SFR) of galaxies are essential 
quests in the study of ga laxy evolution. Whi le optical 
diagnostic diagrams (e.g. [Baldwin et al.l Il981[ ) allows a 
rather clear distinction between star-formation (SF) and 
accretion disk processes, they are limited to - by defini- 
tion - visible components and are at this point extremely 
difficult to apply at high redshifts. The same caveats ap- 
ply to the measurement of SFRs from optical lines. Mid- 
infrared (MIR) spectroscopy offers a potent alternative, 
much less sensitive to interstellar exinction. MIR galaxy 
spectra exhibit an array of features arising essentially 
from (1) a continuous distribution of dust grains, the 
smallest of which (VSGs for Very Small Grains) produce 
the continuum longward of - lO^im (Desert et al. 199(| 
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while larger ones containing silicates produce absorp - 
tion features at 9.7 and 18/im (|Lebofskv fc Riek3ll979( ): 
(2) ionized interstellar gas producing fine-structure lines; 
and (3) molecular gas producing most notably a se- 
ries of broad emission features, most prominent in the 
6 — 17/im range, which were previously referred to as 
"Unidentified Infrared Bands" but are now commonly 
attributed to vibrational emission of l arge polycyclic 
arorn a tic hydocarbon (PAH) ni o lecules (iLeeer fc Pugetl 
[1981 IAllamandola"etaII 11981 iPueet fc Leeed I1989D . 
Rotational lines of m olecular hydrogen are also detected 
(|Roussel et al.l [20071 . and references therein). MIR diag- 
nostics have been devised to unveil the ioniz ing source 
heating these components ( e.g. [V oit 1992b; G enzel et al.l 
119981: [Laurent et all l2000t ISpoon et al.. ,2007i ) and cal- 
ibrations have been proposed to derive SFRs froni 
their luminosities (e.g. 'Ho & Keto"2007'; 'Zhu et al.l|200l 
iRieke et all 120091: [Hernan-Caballero et al. 2009). As 
these calibrations and the resolving power of the various 
diagnostic diagrams vary with galaxy types, it is impor- 
tant to review the MIR spectral properties of well defined 
classes of objects. The Infrared Spectrograph (IRS) on 
board the Spitzer satellite has allowed many such investi- 
gations, building on earlier fu ndamental results from the 
Infrared Space Observat ory (jCesarskv fc Sauvagdl 19991 : 
iGenzel fc Cesarskvll2000f ). Much attention has been de- 
voted to ex treme sources su ch as ULIRGs ([Armus et al.l 
120071: [Fa^ah et al. 2007; D esai et al.l I2007D. starburst 
galax i es dg randl et al. 200 ^, AG"Ns (IWeedman et all 
120051 : iDeo et al., ,200a : .Thompson et all l2009f ) or QSOs 
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Cao et al.ll2008[) IRS observations of the SINGS sample 



Kennicutt et al.|[2003D have also provided many new re- 



sults about the central region of nearby galaxies spanning 
a broad ran ge of physical propert ies ()Dale et al.l 120061 : 
ISmith et al.]l2007t iDale et aljr2009[ ). However few stud- 



ies have yet focused on 'normal' galaxies. Still, questions 
remain open on this seemingly unexciting class of objects. 

Whether VSG or PAH emission can be used to trace 
SF in normal galaxies has been often debated in recent 
vears (.Rousse l et a L 2001; Forster Schreiber et al. ,2004: 
Peeters et aLll2004tlCarzetti et al.ll2007l: iKennicutt et al l 



Resolved observations of star-forming regions 
have shown that the VSG continuum strongly peaks 
inside H ii regions while PAH features dominate 
in photodissociation regions (PDRs) and get weaker 
nearer the core of H ii regions, where the molecules 
are thought t o be destroyed by the intense radia- 
tion fields I'e.g iBoulanger et al.lllQSl: iGiard et al.lll994l: 



Cesarskv et al. 1996; Ve rstraete et al. 119961: iPovich et al.l 
2007; Gordon ct al. 20Qi). However neutral PAH emis- 
sion has recently been reported inside H ii region 
(|Compiegne et al.ll2007[ ). They are also found in the in- 
terstellar medium (ISM), indicating that they must also 
be excited by softer ncar-UV or optic al photons (e.g. 
iLi fc Drainl 12002: Calzctt i et al.l l2007t ). making them 
perha ps better tracer of B stars than of SF (jPeeters et al.l 
[20nl . VSG emission is also observed in the ISM but 
with higher PAH/ VSG surface brightness ratios than in 
SF regions ()Bendo et al.ll200"8 |). Despite much complex- 
ity on small scales however, integrated MIR luminosities 
at 24/j,m and 8/im tracing the VSG and PAH emissions 
respectiv ely, are found t o correlate with Ha luminosi- 
ties (e.g . IZhu et al.l[2008), tho ugh not linearly and with 
scatter ([Kennicutt et al.ll2009l ) leading to uncertain SFR 
estimates. 

An additional source of uncertainty is the common oc- 
curence of AGN in normal galaxies. PAH molecules are 
also thought to get destroyed near the h a rd radiation 
fields of AGNs (iDesert fc Dennefeldl [l98l IWtl fT992al) 
however not totally and as was shown recently from 
IRS spectroscopy, preferentially at short wavelengths 
([Smith et al.ll2007tlO'Dowd et al.ll2009t) . There is in fact 
no a priori reason why PAH emissi on could not be ex - 
cited by UV photons from an AGN ([Farrah et al.ll2007f ). 
This further compromises the use of PAH bands as SFR 
indicators, unless AGNs can be reliably detected in the 
MIR spectra of normal galaxies. 

We have obtained IRS spectra for a sample of 101 nor- 
mal galaxies at z 0.1 with the goal to tackle the above 
issues, making use of additional multi-wavelength (ul- 
traviolet to far-infrared) photometric data and optical 
spectroscopic da ta. The first results o f this survey have 
been reported bv lO'Dowd et al.l ([2009() who analyzed the 
dependence of the relative strength of PAH emission fea- 
tures with optical measures of SF and AGN activity. We 
are pursuing this study by comparing optical and MIR 
diagnostic diagrams to detect AGN presence in normal 
galaxies and by investigating the use of PAH, MIR con- 
tinuum and emission line luminosities as tracer of the to- 
tal IR luminosity and Ha derived SFRs. The sample, IRS 
data and spectral decomposition method are described in 
section 2. Section 3 presents the continuum, PAH and 
emission line properties of the galaxies as a function of 
SF and AGN activity. In particular we analyze the de- 



pendency of PAH equivalent widths with age, metallicity 
and radiation field hardness, as well as the efficiency of 
MIR diagnostics to detect optically classified AGNs in 
these galaxies. We present correlations between the lu- 
minosities of MIR components and the total IR luminos- 
ity in section 4 and between these components and SFR 
estimates in section 5. Our conclusion are summarized in 
section 6. Throughout the paper we assume a fiat ACDM 
cosmology with Hq = 70 km s~^Mpc^^, = 0.3 and 
r^A = 0.7, and a Kroupa IMF (Kroupa .2001.) for SFR 
calibrations. 

2. THE SSGSS SAMPLE 

The Spitzer-SBSS-GALEX Spectroscopic Survey (SS- 
GSS) is a MIR spectroscopic survey of 101 local 
star-for ming galaxie s usin g the Infrared Spectrograph 
(IRS) ([Houck et al.l I2004D aboard the Spitzer Satel- 
lite. The IRS and corollary data a re available at: 



http://www.astro.columbia.edu/ssgss/ 



2.1. The parent sample 

The sample is drawn from the Lockman Hole region 
which has been extensively surveyed at multiple wave- 
lengths. In particular UV photometry from GALEX 
(1500 and 2300A), optical imaging and spectroscopic ob- 
servations from SDSS and infrared photometry (IRAC 
and MIPS channels) from Spitzer (SWIRE) are available 
for all SSGSS galaxies. The redshifts span 0.03 < z < 
0.21 with a mean of 0.09 similar to that of the fuU SDSS 
spectroscopic sample. The sample has a surface bright- 
ness limit of 0.75 MJy sr"-'^ at 5. 8/im and a flux limit 
of 1.5mJy at 24/im. Due to these cuts the sample does 
not contain very low mass/low metallicity /low extinction 
galaxies (9.3 < M© < 11.3, 8.7 < log(0/H) + 12 < 9.2 
and 0.4 < A^a < 2.3) but it was selected to cover the 
range of physical properties of 'normal' galaxies. 

These galaxies are divided up into 3 categories: star- 
formation dominated galaxies (referred to as 'SF galax- 
ies'), composite galaxies (SF galaxies with an AGN com- 
ponent) and AGN dominated gala xies, according t o thei r 
location on the "BPT diagram" ([Baldwin et aLill98lD . 
which shows [N ii]A6583/Ha (a proxy for gas phase 
metallicity) against [O iii]A5007/H/3 (a measure of the 
hardness of the radiation field). AGNs (all Seyfert 2's 
in our sample) are isolated by the theoretical bound- 
ary of iKewlev et al.l ((200 1() while SF galaxies and com- 
po site galaxies are s e parat ed by the empirical boundary 
of lKauffmann et all l|200l . In all following figures, SF 
galaxies are represented as black dots, composite galaxies 
as pink stars and AGNs as open red triangles. We note 
that this optical classification may miss obscured AGNs 
contributing to the MIR emission. 

The location of the sample in the NUV — r color versus 
r-band absolute magnitude diagram is shown in Figure 
d] (left panel) with the volu me density contou rs of the 
underlying local population ([Wvder et al.! '2007V Galax- 
ies in this diagram separate into two well-defined blue 
and red sequences that become redder with increasing 
luminosity. The red sequence tend to be dominated by 
high surface brightness, early type galaxies with low ra- 
tios of current to past averaged SF, while the blue se- 
quence is populated by morphologically late-type galax- 
ies w ith lower surface brig htness and on-going SF activity 
(e.g. IStrateva et aLll2001f) . The color variation along the 
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Fig. 1. — Left: NUV — r color vs r-band magnit ude diagram showi ng the location of SSGSS galaxies with respect to the underlying 
local population shown as volume density contours l lWvder et al.ll2b07l) . Galaxies separate into two well-defined blue and red sequences. 
In this and all following figures, star-forming galaxies are represented as black dots, composite galaxies as pink stars and AGNs as open 
red triangles. The SSGSS sample is dominated by blue sequence galaxies with a small fraction (mostly AGNs) on the red sequence. Right: 
Infrared color-color diagram: fs/f24 versus /70//16O flux ratios. 'Normal' quiescent SF galaxies are found towards the bottom right corner 
(stronger PAH emission and cooler dust in the FIR) while starburst galaxies populate the top left corner (strong hot dust continuum in 
the MIR and warm dust emission in the FIR) . The Spitzer data for NGC 35321, NGC 337 and Mrk 33 are taken from lDale etaLl 1 I2OOI) . 



blue sequence is due to a combination of dust, SF his- 
tory and metallicity (jWvder et al.ll2007|) . Unsurprisingly 
given its selection criteria, our sample is dominated by 
blue sequence galaxies, although a small fraction (mostly 
AGNs) are found on the red sequence. 

The right panel of Figure [1] shows the distribution of 
the sample in the fs/ f2i - /70//16O plane, where fs/ fiA 
is the 8/im to 24/im restframe flux ratio (4*'' IRAC band 
to MIPS band) and /70//16O is the TO^im to 160^m 
restframe flux ratio (2"'^ to S'''^ MIPS bands). The in- 
frared fc-corrections are described in Section 4. This fig- 
ure illust rates the range of IR p roperties covered by the 
sample ( da Cunha et al.l ()2008[ ) and references therein). 
Starburst galaxies tend to populate the top left corner 
(strong hot dust continuum in the MIR and warm dust 
emission in the FIR) while more quiescent SF galaxies 
move down towards the bottom right corner (stronger 
PAH emission and cooler dust in the FIR). The Spitzer 
data for NGC 3 5321, NGC 337 and Mrk 33 are from 
IDale et al.l ([2007I) . 

2.2. The IRS Spectra 

Low resolution spectroscopic observations were ac- 
quired for the full sample using the Short-Low (SL) and 
Long Low (LL) IRS modules, covering 5/im to 38/im 
with resolving power ^ 60 to 127. High resolution spec- 
tra were obtained for the 33 brightest galaxies using the 
Short-High (SH) IRS module, covering lO^m to 19.6^m 
with a resolving power of ^ 600. A detailed descrip- 
tion of the data acquisition and reduction can be found 
in O'Dowd et al. (in preparation). In short, standard 
IRS calibrations were performed by the Spitzer Pipeline 
version S15.3.0 (ramp fitting, dark substraction, droop, 
linearity correction, distortion correction, flat fielding, 
masking and interpolation, and wavelength calibration). 




Fig. 2. — Two example spectra with best fits {ulv on the 
y-axis in uriits of 10^ ^Jy Hz). The fits are outputs of PAHFIT 
llSmith et al.|[2007^ . The red lines fit the continuum, the purple 
lines fit the PAH features and the dotted green lines fit the emis- 
sion lines. The main PAH features are indicated in purple and the 
main emissions lines in green. The blue curves show - from left 
to right — the filter responses of the IRAC bands at 6 and 8/im, 
of the IRS blue Peak Up band at I6fj,m and of the MIPS band at 
24/jm bands. The galaxy in the top panel is a SF galaxy with no 
silicate absorption {tx = 0); the galaxy in the bottom panel is an 
AGN with strong silicate absorption features at 9.7fim and ISfim 
(tx = 2.5). The extinction (1 — e~'^-^)/T>, is shown as the dotted 
line in arbitrary units. 



Sky subtraction was performed manually with sky frames 
constructed from the 2-D data frames, utilizing the shift 
in galaxy spectrum position between orders to obtain 
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clean sky regions. IRSCLEAN (vl.9) was used to clean 
bad and rogue pixels. SPICE was used to extract 1- 
D spectra, which were combined and stitched manually 
by weighted mean. After rejecting problematic data, 
the final sample consists of 82 galaxies (56 SF galaxies, 
19 composite galaxies and 7 AGNs) with low resolution 
spectra, of which 31 (23 SF galaxies, 6 composite galaxies 
and 2 AGNs) have high resolution spectra as well. 

2.3. PAHFIT decomposition 

We us e the PAHFIT sp ectral decomposition code 
(vl.2) of lSmith et all pOOT ) (hereafter S07) to fit each 
spectrum as a sum of dust attenuated starlight contin- 
uum, thermal dust continuum, PAH features and emis- 
sion lines. The absorbing dust is assumed to be uniformly 
mixed with the emitting material. The code performs a 
fitting of the emergent fiux as the sum of the following 
components (Eq. 1 in S07): 



M 



^(A/Ao)2 



(1 







(1) 

where B^, is the backbody function, = 5000 K 
is the temperature of the stellar continuum, Tm — 
35, 40, 50, 65, 90, 135, 200, 300 K are 8 thermal dust con- 
tinuum temperatures, the Irii^) consist of 25 PAH emis- 
sion features modeled as Drude profiles and 18 unresolved 
emission lines modeled as Gaussian profiles, and t\ is the 
dust opacity, normalized at Ao — 9.7/im. The specifics of 
these components are described in S07. The Drude pro- 
file, which has more power in the extended wings than a 
Gaussian, is the theoretical profile for a classical damped 
harmonic oscillator and is thus a natural choice to model 
PAH emission. Some of the PAH features are modeled 
by several blended subfeatures, most prominently the 
PAH complex at 7.7fim which is modeled by a combina- 
tion of 3 Drude profiles and the PAH complex at 17/im 
modeled by 4 such profiles. The continuum components 
have little significance individually, it is their combina- 
tion that is meant to produce a physically realistic con- 
tinuum. We find that the stellar continuum is negligible 
for most galaxies, which is probably not surprising since 
it is practicall y unconstrained. The theoretical value of 
0.232 X /s.e^m ()Helou et al.ll2004[) for the stellar contribu- 
tion to the 8/im band is ~ 20%, however it is probably an 
upper limit since the 3.6/im flux may be contaminated by 
the 3.3/im PAH feature for galaxies at z ~ 0.1. We also 
find that silicate absorption is negligible (rg.y < 0.1) for 
63% of the sample, however 3 out of 7 AGNs are among 
the galaxies showing the strongest absorption features. 
From the best fit decompositions^ we compute the fluxes 
of the PAH features, emission lines and continuum at 
various points corrected for silicate absorption as well 
as the total restframe and observed fluxes in the MIR 
Spitzer bands. The fluxes of the main MIR components 
used in this paper are listed in Table [T] 

Figure [5] shows two examples of our IRS spectra with 
best fit decomposition from PAHFIT (restframe). The 
galaxy in the top panel is a typical SF galaxy (#30); 
that in the bottom panel is an AGN (#93) showing the 
strongest silicate absorption features at 9.7/im and 18/im 

^ corrected for the PAHFIT vl.2 (1 + z) overestimate. 
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Fig. 3. — Observed PAH fractions in the 8/im IRAC band and 
16/jm IRS band as a function of redshift. Symbols are as described 
in Figure [T] The 8/im flux is dominated by PAH emission for most 
galaxies while the 16/tm flux is dominated by continuum emission 
except for the two highest redshift sources. 



in the sample (rg.7 — 2.5). The extinction (1 — e^'^^)/T\ 
is shown as the dotted line in arbitrary units. 

Figure [3] shows the observed fractions of PAH emission 
in the 8/im IRAC band (often used as a proxy for PAH 
emission) and in the 16/im IRS band as a function of red- 
shift. For our local sample the 8/im IRAC channel picks 
up the 7 .IpLia. PAH complex, plus the 6.2/im PAH feature 
for galaxies ai z > 0.05 and the 8.6/xm PAH feature for 
galaxies at z < 0.05. Both the observed and restframe 
fiuxes in this band are largely dominated by PAH emis- 
sion for most SF and composite galaxies. The continuum 
dominates only for 1 AGN. The 16/im IRS Peak-Up band 
collects photons from the 17//m PAH complex plus other 
smaller PAH features around 14/im and the large 12.7/im 
complex for galaxies at 2 > 0.06. The observed 16/im 
fiux includes more PAH emission than the restframe flux, 
which is continuum dominated (~ 70%) for all galaxies. 
The observed 24/xm MIPS channel is vastly dominated 
by the continuum for all sources. The highest PAH con- 
tribution (18%) comes from the redshifted 18.92/im PAH 
feature and the red wing of the 17/im PAH feature for the 
highest redshift object {z = 0.217). PAH emission starts 
to dominate the 24/im channel for galaxies at z > 1. 

2.4. Aperture corrections 

The SL and LL IRS modules have slit widths of 3.6" 
and 10.5" respectively, while the mean angular size of the 
sample is 10" . The corrections applied to stitch the two 
modules together in the overlap region (14.0 — 14.5/im) 
are explained in detail by O'Dowd et al. (in prepara- 
tion). Wavelength dependent aperture effects also arise 
from the wavelength dependent PSF (increased sampling 
of the central regions of extended galaxies with increas- 
ing wavelength). To remedy these effects, we compute 
spectral magnitudes in the MIR Spitzer bands both from 
the data and from the PAHFIT SEDs using the Spitzer 
Synthetic Photometry cookbook. We flnd excellent agree- 
ment between data and fits except in the 6/im IRAC 
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Fig. 4. — Aperture corrections in magnitude at 6, 8, 16 and 24/im 
as a function of r-band Petrosian diameter in arcseconds. The 
dotted lines mark the average aperture corrections. The vertical 
lines show the slit widths of the SL and LL modules (upper and 
lower panels respectively). 



bands, the data being noisy and the fits unreliable below 
5.8/im. For this reason we do not make use of fluxes in 
this part of the spectrum. The difference between the 
PAHFIT spectral magnitudes and the photometric mag- 
nitudes are used as aperture corrections at the effective 
wavelengths. These corrections are shown in Figure |4] 
as a function of r-band petrosian diameter and listed in 
Table [2] The vertical lines in the upper and lower panels 
show the slit widths of the SL and LL modules respec- 
tively. It is clear that flux is lost at 8/im, however at 
longer wavelengths we do not find that much flux is lost 
even when the optical petrosian diameter is larger than 
the slit width, which we attribute to the larger PSF. 
Corrections at intermediate wavelengths are obtained by 
interpolation. The mean corrections are ~ 1.2 mag at 
6/im, ~ 0.5 mag at 8/im and < 0.1 at 16 and 24//m. In 
the following, all MIR luminosities computed from the 
PAHFIT decomposition (PAH, continuum and emission 
line luminosities as well as total rcstframc luminosities in 
the Spitzer bands) are corrected for aperture as described 
in this section. 

3. MIR SPECTRAL PROPERTIES 

Figure [5] shows the mean spectra of our SF galaxies 
(solid line), composite galaxies (dotted line) and AGNs 
(d ashed line) as well a s the average starburst spectrum 
of iBrandl et al.l ()2006l ) (dot-dash), normalized at 10/im. 
The transition from starburst to SF galaxy to AGN 
is associated with a declining continuum slope, most 
dramatic between the starburst spectrum and the nor- 
mal SF spectrum. Indeed H ii regions and starburst 
galaxies are found to exhibit a steep risi ng VSG contin 



uum component longward of ^ 9p,m ( e.g.'^Cesarsk v et al 
1996HLaurent et al.ll2000HDale et all 2001: Pcctcr s et al 



2004^ . The transition is also marked by decreased 



Fig. 5. — The mean spectra of SF galaxies (solid line), com- 
posite galaxies (dotted line) and AGNs (dashed line) normalized 
at lO^m. The dot-dashed spectrum is the average starburst spec- 
trum of Brandl et al. (2006). The transition from starburst to SF 
galaxy to AGN is marked by a declining continuum slope, decreased 
[Ne n]12.8Atm and [S ni]18.7/im and enhanced [O lv]25.9Atm. The 
AGN and starburst spectra also show depleted PAH emission at 
low wavelength compared to the SF spectrum. 




[Ne ii] 12. 8/im and [S iii]18.7/im line emis sion and en- 
hanc ed [O iv]25.9/im line emission (e.g. iGenzel et alJ 
I1998D . The AGN spectrum, and to a lesser extent the 



X [,um] 

Fig. 6. — Top: The mean PAH component of 'young' SF galaxies 
with 1.1 < Dn('iOOO) < 1.3 (solid line) and 'old' SF galaxies with 
1.3 < D„(4000) < 1.6 (dashed Une). Bottom: The mean PAH 
component of SF galaxies (solid line) and AGNs (dashed line) in the 
1.3 < D„(4000) < 1.6 range where both types have similar mean 
Dn(4000) ~ 1.4. The spectral components are normalized by the 
peak intensity of the 7.7/xm feature. The main difference between 
the two pairs is enhanced PAH emission at large wavelengths with 
respect to the 7.7^m feature. 



starburst spectrum, show weaker PAH emission at low 
wavelength than the SF spectrum, an effect attributed 
to the destruction of PAHs in intense far-UV radiation 
fields. 

3.1. PAH features, continuum and emission lines 
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Fig. 7. — MIR color indices a(8, 16) and o{16, 24) as a function 
of D„(4000), [O lll]A5007/H/3 ratio and /70//16O restframe colors. 
Symbols are as described in Figure [T] Although older galaxies and 
AGNs have shallower continuum slopes on average, little corre- 
lation is found with Dn(4000) nor radiation field hardness. The 
correlation with FIR color for SF galaxies may reflect a sequence 
in the peak wavelength of the dust SED: the MIR slope steepens 
at it gets closer to the peak while /70 //16O increases. 



PAHFIT allows us to compare the different spectral 
components of different galaxy types separately. The 
top panel of Figure [S] shows the average PAH compo- 
nent of 'young' SF galaxies with 1.1 < D„(4000) < 1.3 
(< D„(4000) >= 1.2) and that of 'old' SF galaxies with 
1.3 < L'„(4000) < 1.6 (< D„(4000) >= 1.4). The 4000A 
break D„(4000) (|Balogh et al.lll998[ ) is a measure of the 
average age of the stellar populations. The separating 
value is simply the median of the distribution. The bot- 
tom panel shows the average PAH components of SF 
galaxies and AGNs in the 1.3 < £)„(4000) < 1.6 range 
where both types have similar mean D„(4000) ^ 1.4 
(there are only 4 AGNs in that bin) . All spectra are nor- 
malized by the peak intensity of the 7.7//m feature. The 
main difference between the pairs in both panels is en- 
hanced PAH emission at large wavelengths with respect 
to the 7.7^m feature, i.e. an increase in the ratio of high 
to low wavelength PAHs associated with both AGN pres- 
ence and increased stellar population age. This increase 
is most pronounced in the lower panel (AGN versus SF) 
where a decrease in the 6.2/xm feature with respect to 
the 7.7/Km feature is also noticeable. The variations in 
PA H ratios in this sample have been thoroughly studied 
by lO'Dowd et al.l (2009) and shown to be statistically 
significant. These variations are much more dramatic 
for AGNs with harder radiation fields than those in the 
present sample (e.g. S07, their Figure 14). They can be 
attributed to a change in the fraction of neutral to ion- 
ized PAHs responsible for the high and low wavelength 
features respectively, and/or to the destruction by hard 
radiation fields in AGNs of the smallest PAH grains emit- 
ting at low wavelengths (S07, and references therein). 
The variations of PAH strengths with age, metallicity 
and radiation field hardness are explored in more detail 
in the next section. 

We define the continuum slope or MIR color index be- 
tween wavelengths Ai and A2 as: 



a(Ai, A2) = 



log [/r'*(A2)//r"*(Ai)] 



(2) 



[Nell] 
iJh [Nelll] 



o 



51^ 



log(A2/Ai) 

where /™"* (A) is the continuum component of Eq. [T] at 



I [Jy*Hz] HIGH RES 

Fig. 8. — Comparison between the low and high resolution line 
fluxes of [Ne ll]12.8/im (in black) and [Ne in]15.5^m (in blue) for 
a subsample of 31 galaxies. Excluding the objects marked with a 
red cross, the rms of the correlation is 0.22. 



A corrected for silicate absorption. This would be the 
index /3 of a continuum spectrum of the form I^, oc A*^. 
Figure [7] shows a(8, 16) and a(16, 24) as a function of 
D„(4000), [O iii]A5007/H^ and the restframe /70//16O 
color (see Section 4 for details on the fc-corrections) . As 
discussed above and shown in Figure [SJ the mean MIR 
slope is found to steepen fro m quiescent galax ies to star- 
burts of increasing activity (Dale ct al."2001") and to be 
shallower for AGNs (e.g. Gcnzcl & Ccsarsky 2000). How- 
ever our indices span a significant range (~3 dex) with 
little correlation with the age of the stellar populations or 
radiation field hardness. Older galaxies (D„ (4000) > 1.6) 
do tend to populate the low end of the distribution (i.e. 
have shallower slopes) in both cases, as do AGNs in the 
red part of the spectrum however a flatter continuum 
could not be used as a criterion to sepa rate AGNs from 
SF ga laxies, as previously reported by IWeedman et al.l 
()2005l ). The correlation with FIR color for SF galaxies 
is more striking, especially at longer MIR wavelengths. 
This may be expected if the peak of the dust SED (a 
blackbody modified by the emissivity) is located short- 
ward of ~ 100/im. In this case as the peak wavelength 
decreases, the MIR continuum slope gets closer to the 
peak and therefore steepens while /70//16O increases. 

Finally we look at variations in the emission line com- 
ponents. The lines modeled by PAHFIT in the low res- 
olution spectra are meant to provide a realistic decom- 
position of the blended PAH features and the contin- 
uum (S07) but the spectral resolution is of the same 
order as the FWHM of the lines. Figure [8] shows the 
comparison between the high and low resolution fluxes 
of the [Ne ii]12.8/.tm and [Ne iii]15.5/Ltm lines (black 
and blue error bars respectively) for the subsample ob- 
served with the SH module. The high resolution lines 
were also measured using PAHFIT with the default set- 
tings. We make no attempt at aperture correction on 
this plot. Excluding the 3 extreme error bars among the 
[Ne iii]15.5/xm fluxes at high resolution and the outlier 
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Fig. 9. — Top: The difference - A/^ - between the mean emission 
line component of 'young' SF galaxies with 1.1 < D,i(4000) < 1.3 
and that of 'old' SF galaxies with 1.3 < D„(4000) < 1.6. Bottom: 
The difference between the mean emission line component of SF 
galaxies and that of AGNs in the range 1.3 < D„(4000) < 1.6. 
The spectral components are normalized by the total flux in the 
16/im IRS band. The most significant features are decreased 
[Ne lll]15.5/im and [S ni]18.7/im and increased H25(0) in older SF 
galaxies with respect to younger ones, and increased [O iv]25.9^m 
and H2S(1) in AGNs along with diminished [Ne n]12.8/jm and 
[S in]18.7/im with respect to SF galaxies. 



among the [Ne ii]12.8/im fluxes (marked as red crosses in 
Figure H]), the fitting procedure at low resolution recovers 
the high resolution fiuxes with an rms of 0.22 dex, a rea- 
sonable estimate considering the factor of 10 difference 
in spectral resolution. In particular the PAH contami- 
nation for the [Ne ii]12.8/im line does not seem to be a 
significant problem in the SL data using PAHFIT. For 
the purpose of the present statistical analysis we use the 
low resolution line measurements which are available for 
the full sample and over the full range of wavelengths. 
We refer to O'Dowd et al. (in preparation) for a detailed 
comparison between the high and low resolution data. 

The top panel of Figure [3] shows the difference - A/^ - 
between the average emission line component of 'young' 
SF galaxies (< Ai(4000) >= 1.2) and that of 'old' 
SF galaxies (< i:)„(4000) >= 1.4) as defined earlier, 
while the bottom panel shows the difference between 
the mean emission line component of SF galaxies and 
that of AGNs in their overlapping range of D„(4000) 
(1.3 < i:>„(4000) < 1.6). The spectral components 
were normalized to the total fiux in the 16/im IRS band. 
Among the most significant features are the decreased 
[Ne iii]15.5/xm and [S iii]18.7/im lines and increased 
112^(0) line in older SF galaxies with respect to younger 
ones, and the strong increase in [O iv]25.9/im line emis- 
sion in AGNs along with diminished [Ne ii]12.8/im and 
[S iii]18.7/im emission with respect to SF galaxies. The 
Il2S(l) molecular line is also enhanced in AGNs. A 
strong excess of H2 in many Seyferts and LINERS has 
been reported by Rous sel et al.l (|2007[ ). suggesting a dif- 
ferent excitation mechanism in these galaxies. H2 line 
emission is studied in more detail in Section 5.4. 

While low excitation lines such as [Ne ii]12.8/im and 
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Fig. 10.— Top: The [Ne lll]15.5Afm/H2S'(0) line ratios as a 
function of D„(4000). Bottom: The [Ne ll] 12.8/im/[0 iv]25.9/im 
line ratios as a function of [O ill] A5007/H/3. Symbols are as de- 
scribed in Figure [T] Ratios with extremely large errors are not 
shown. The pearson coefficients of the correlations are indicated 
in each panel for the full sample. The most significant correlation 
is found between [Ne ll] 12.8/im/[0 iv]25.9/im and [O III] A5007/H/3 
for the subsample of composite galaxies and AGNs (r = —0.80). 



[Ne iii]15.5/xm can be excited by hot stars as well as 
AGNs (they are detected in all but 1 spectrum for 
[Ne ii]12.8/xm, all but 3 spectra for [Ne iii]15.5^m), 
the high excitation potential of the [OIV]25.89/zm line 
(54.9eV) (the brightest such line with [NeV]14.21^m 
in the MIR ) usua l ly links it to AGN activity (e.g. 
Genzel et all Il998t ISturm et al] 120021: iMelendez etall 



20081) . However it has also been attributed to star- 



burst related n iechanisms (jSchaerer fc Stasihskal 119991 : 
iLutz et anil998l ) and indeed detected in starburst galax- 
ies or regions (|Lut z et al. 1998; Bcirao et al. 200^ 
lAlonso-Herrero et al.l l2009f ). It is detected in 73% of 
our 'pure' star-forming galaxies (63% of the composite 
galaxies) while undetected in 1 out of 7 AGNs. It may 
also be that three quarters of our SF galaxies harbor 
an obscured AGN not detected in the optical. While 
the non detection of [OIV]25.8 9Atm in AGNs has also 
been known to happen (e.g. iWeedman et al.l [20051 ) . the 
one AGN spectrum in our sample without [OIV]25.89/im 
(#63) is particularly noisy and the presence of the line, 
even significant, cannot be ruled out. 

The top panel of Figure [10] shows the 
[Ne iii]15.5^m/H2S'(0) ratios as a function of D„(4000). 
The Pearson coefficient of the correlation is indi- 
cated in the top right corner. The trend is mild, 
and milder still for the [S in]18.7/zm/H2S'(0) ratios. 
Much more significant is the correlation between 
[Ne ii]12.8/.tm/[0 iv]25.9fim and [O iii]A5007/H/3 
shown in the bottom panel. The correlation for 
[S iii]18.7/im/[0 iv]25.9/im is somewhat less significant 
but both ratios notably decrease with increasing radi- 
ation field hardness for composite galaxies and AGNs 
(the Pearson coefficient for this subsample is r = —0.80). 
Ratios of high to low excitation emission lines have 
long been used to characterize the dominant source of 
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Fig. 11. — The [Ne lll]15.5/xm/[Nc ll]12.8/im emission line ratios 
of star-forming galaxies as a function o f metallicity. Th e inset 
shows the low met allicity d ata points of lO'Halloran et al.l II2006I ) 
(open squares) and IWu et a l. (2006) (green error bars and lower 
limits), with our dynamic range shown as the dotted box in the 
bottom right corner. 
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Fig. 12. — The 6.2/im PAH equivalent widths (EWs) computed 
using Eq. [3]with the PA HFIT decomposition paramete rs compared 
to the EWs computed bv lSargsvan &: Weedmanl II2009I) assuming a 
single gaussian on a linear continuum between 5.5/jm and 6.9/xm. 
Symbols are as described in Figure [T] 



ionization in galaxies (e.g. IGenzel et allll998[ ). We come 
back to this point in Section 3.3. 

The [Ne ili]15.5/zm/[Ne il]12.8yLtm line ratio is also ex- 
pected to be sensitive to the hardness of the radiation 
field, however we find no correlation between this ratio 
and [O iii]A5007/H/3 in our sample. We do find a trend 
with metallicity despite the very narrow metallicity range 
of our sample, as shown in Figure fTD for the SF subsam- 
ple. Indeed [Ne n]12.8fim has been shown to be the dom- 
inant ionization species in H ii region at high metallicity 
while [Ne iii]15.5/xm takes over in regions of lower den- 
sity and higher excitation such as low mass, low metal- 
licity galaxies (jO'Halloran et al.ll2006t IWu et al.l [20061 ) . 
The inset shows a larger scale v ersion of this figure with 
low metallicity data points from lO'Halloran et all (|2006f ) 
(open squares) and IWu et al.l ()2006[) (green error bars 
and lower limits). Our dynamic range is represented as 
the dotted box in the bottom right corner. 

3.2. PAH Equivalent Widths 

We compute equivalent widths (EWs) as the integrated 
intensity of the Drude profile(s) fitting a particular PAH 
feature, divided by the continuum intensity below the 
peak of that feature. Using Eq. 3 from S07 for the 
integrated intensity of a Drude profile, the EW of a PAH 
feature with central wavelength A^, full width at half- 
maximum FWHMr (as listed in S07, their Table 3), and 
central intensity br (PAHFIT output), can be written as: 

where /™"*(Ar) is the continuum component of Eq. [TJ 
This definition is different from that of S07 in PAH- 
FIT which computes the integral / (/^^^//^°"*)dA in 
the range A^ ± 6 x FWHMr- In the case of the 7.7 mi- 
cron feature whose FWHM is large and extends the limit 
of the integral to regions beyond the IRS range where the 



continuum vanishes arbitrarily, the profile weighted aver- 
age continuum is used. Despite this caveat, both meth- 
ods agree within 10% and the discrespancies virtually 
disappear when increasing the limits of the integral for 
all other PAHs^". However the EWs measured as above 
differ significantly from those estimated with the spline 
method, which consists in fitting a spline function to the 
continuum from anchor points around the PAH feature, 
and a Gaussian profile to the continuum-subtracted fea- 
ture. This method yields considerably smaller EW values 
as it assigns a non negligible fraction of the PAH flux ex- 
tracted by PAHFIT to the continuum. Figure |T2] shows 
our EWs (Eq. H against the 6 . 2Mm PAH EWs com- 
puted by iSargsvan fc Weedmanl ([200l for the SSGSS 
sample assuming a single gaussian on a linear continuum 
between 5.5/im and 6.9/im. Their published sample is 
restricted to SF galaxies defined as having EW(6.2/.tm) 
> 0.4/xm (jWeedman fc Houc"3l2009l ). The measurements 
for the remaining galaxies were kindly provided by L. 
Sargsyan. Their formal uncertainty is estimated to be 
~ 10%. The two methods are obviously strongly diver- 
gent. The spline EWs strongly peak around a value of 
0.6/im with no apparent correlation with the PAH- 
FIT estimates, which reach ^ 15 nm and can be up to 25 
times larger than the ISargsvan fc Weedmanl (|2009l) val- 
ues. Our EWs for the main PAH features are listed in 
Table O 

The strength of a PAH feature depends on several 
interwined properties of the ISM: metallicity, radia- 
tion field hardness, dust column density, size and ion- 

In the process of making these comparisons, we discovered 
two bugs in PAHFIT: 1/ the code was mistakenly calling gaussian 
profiles instead of Drude profiles to compute the EW integral, thus 
underestimating EWs by ~ 1.4, and 2/ it was applying silicate 
extinction to the continuum while using the extinction corrected 
PAH features (according to Eq. 1, both components are equally 
affected by the extinction term). These bugs are being corrected 
(J.D. Smith, private communication). 
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Fig. 13. — The equivalent widths of the main PAH features at short wavelengths (6.2, 7.7 and 8.6/xm, left panel) and at long wavelengths 
(11.3, 12.7 and 17/.tm, right panel) as a function of [O III] A5007/H/3. Symbols are as described in Figure [T] The Pearson correlation 
coefficients r are indicated at the top right of each panel. The dotted lines in the left panel show the limits that best isolate SF galaxies 
(EW(6.2/im) > 1/xm, EW(7.7/.im) > 4/im and EW(8.6/im) > Ifim). The dotted lines in the right panel are approximate lower limits for the 
SF population EW(11.3/^m) > l.S^m, EW(12.7/im) > 0.9/im and EW(17/im) > O.S^im). AGN EWs become increasingly undistinguishable 
from those of SF galaxies towards longer wavelengths. 



ization state distributions of the dust grains (Dale et 
al. 2006 and references therein). In particular it is 
shown to be reduced in extreme far-UV radiation fields, 
such as AGN-dom inated en vironment (Gcnzc l et alJ 
Il998t ISturm et all 2000; Wcc dman et al . 2005) , near 
the si tes of SF (iGeballe et alTTl989l: ICesarskv et alJ 



1996t iTacconi-Garman et alJ l2005t iBeirao et all 120061: 



Povich et all 120071: iGordon et al.' '2008*) or in very low 
metallicity environmen ts (Dwck 2005: Wu et al.l 120051 : 
lO'Halloran et alJ 120061: iMadden et al.1 [20061), where the 
PAH molecules are thought to get destroyed (e.g. IVoitI 
[i992a). 

Figure [13] shows the EWs of the main PAH features 
as a function of [O iii]A5007/H/3. The Pearson correla- 
tion coefficients r are indicated at the top right of each 
panel. AGNs do exhibit noticeably smaller EWs than SF 
galaxies at short wavelengths (6.2, 7.7 and 8.6//m, left 
panel), however seemingly uncorrelated with radiation 
field hardness. The range of EWs spanned by AGNs be- 
comes increasingly similar to that of SF galaxies towards 
longer wavelengths (11.3, 12.7 and 17fim, right panel) 
while at the same time a correlation seems to appear 
with radiation field hardness. The Pearson coefficients 
for the AGN population alone at long wavelengths are 
-0.97, -0.89 and -0.81 respectively from top to bottom, 
though admittedly they are boosted by the rightmost 
data point. A larger sample of AGNs is needed to con- 
firm this correlation. PAH strength remains largely in- 
dependent of radiation field hardness for SF and com- 
posite galaxie s. The se results complement the anal- 
ysis of ,0'Dowd et al.l ((2009) who found a correlation 
between the long-to-short wavelength PAH ratios and 
[O iil]A5007/H/3 in AGNs. These trends are consistent 
with the selective destruction of PAH molecules in the 
hard radiation fields of these sources ([O iii]A5007/H/3 
> 1.5). The EW trends or lack thereof in Figure [13] 



suggest that the smallest PAH molecules effective at 
producing the short-wavelength PAH features get de- 
stroyed first, near an AGN, while the larger molecules 
producing the larger wavelength PAHs require increas- 
ingly harder radiation fields for th eir PAH strength to 
drop below that of SF galaxies. iDesert fc Dennefeldl 
(|1988D first suggested that the absence of PAHs could 
be taken as evidence for the presence of an AGN. Weak 
PAH emission has since often been used to discriminate 
between photoionization and accretion disk processes. 
However the common boundaries for a 'pure starburt', 
e.g. EW f7.7um)> 1 (iLutz et al.lll998D or EW(6.2/.tm) 
> 0.4/xm (jWeedman fc Houckll2009i r are significantly too 
weak here, due to the different method we use to com- 
pute the equivalent widths as shown above. Based on 
the PAHFIT decomposition, SF galaxies would be best 
isolated by EW(6.2^m) > l/im, EW(7.7^m) > 4^m or 
EW(8.6/-tm)>l/im, the latter two criteria being more ac- 
curately determined in our sample. Those limits are 
shown as dotted lines in the left panel of Figure [T51 
The two SF exceptions below the 7.7 fim and 8.6yLtm EW 
limits (#32 and #74) happen to have very strong sili- 
cate absorption parameters (rg.y =1.8 and 2.33) and still 
very distorted absorption-corrected continua compared 
to the rest of the sample. The dotted lines in the right 
panel are approximate lower limits for the SF popula- 
tion (EW(11.3/xm) > 1.8/zm, EW(12.7/xm) > 0.9^*™ and 
EW(17^m) > O.Sfim). It is clear that the AGN popu- 
lation becomes increasingly difficult to isolate based on 
EW alone in the red part of the spectrum. 

Figures [T4l and [T5] show the EWs of the main PAH fea- 
tures as a function D„(4000) and [N li]A6583/Ha respec- 
tively. The EWs at short wavelengths show a mild down- 
ward trend with increasing age (or decreasing SF activi- 
ties) while they become independent of it at long wave- 
lengths. This again is consistent with the correlations 
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Fig. 14. — The equivalent widths of the main PAH featur es a t short and long wavelengths (left and right panels respectively) as a function 
of D„(4000). Symbols and lines are as described in Figure [13] 
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Fig. 15. — The equivalent widths of the main PAH features at short and long wavelengths (left and right panels respectively) as a 
function of [N ll]A6583/Ha. Symbols and lines are as described in Figure [13] The short wavelength PAH EWs significantly decrease with 
[N ii]A6583/HQf. 



between the long-to-short wavelength P AH ratios and 
Pn(4000) or Ha equivalent width found bv lO'Dowd et alj 
(|2009i ). The short wavelength EWs decrease more no- 
tably with increasing [N ii]A6583/Ha ratios, which of 
course are related to D„ (4000) but appear to be the prop- 
erty that most uniformally and significantly affects the 
sample as a whole. Metallicity and SF activity are known 
to affect PAH strength, however, as mentionned earlier, 
previous studies have demonstrated the opposite effect, 
namely a decrease in PAH strength at very low metal- 
licity and in intense SF environment. These trends thus 
make normal blue sequence galaxies the sites of maxi- 
mum PAH strength. 
Other than PAH destruction, another cause of decreas- 



ing PAH strength at low wavelengths may be a stronger 
continuum whose strength may depend on the above pa- 
rameters. A short wavelength continuum (3-10/im) has 
been observed in AGNs, which is attributed to very hot 
dust heated by their intense radiation fields (Laurent et 
al. 2000 and references therein), however the continuum 
slopes of AGNs in our sample largely overlap those of 
the SF population (Figure [T]). Figure [TC] shows the EW 
of the 7.7/im feature as a functio n of continuurn . slope , 
a diagnostic diagram proposed bv lLaurent et al.l ()2000l ) 
to distinguish AGNs from PDRs and H ii regions. A 
clear trend is seen for SF and composite galaxies, sug- 
gesting that decreased PAH strength in normal SF galax- 
ies may be at least partly due to an increased continuum 
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Fig. 16. — The equivalent width of the l.T/im feature as a 
function of contin uum slope (a diagnostic diagram proposed by 
ILaurent et al.l 11200 0)). Symbols are as described in Figure [T] A 
clear trend is seen for SF and composite galaxies while AGNs ap- 
pear more randomly distributed. AGNs are expected to populate 
the lower left side of the plot, PDRs the lower right corner and H II 
regions the upper left corner of the diagram (.Laurent et al...2000.) . 



at low wavelength, which is itself loosely correlated with 
D„(4000) (FigureEl). However the EWs of AGNs appear 
quite independent of their continuum slope, supporting 
the PAH destruction scenario. In this diagram, AGNs 
are expected to populate the lower left side of the plot 
(shallow slopes and weak PAH features) , PDRs the lower 
right corner (shallow slopes and strong PAH features) 
and H ii regions the upper left corner of the diagram 
(steep slopes and weak PAH features) . Our SF sequence 
is qualitatively similar to the location of quiescent SF 
regions on the Laurent et al. diagram (their Figure 6), 
which are modeled by a mix of PDR and H ii region 
spectra, plus an AGN component towards the lower left 
corner where composite galaxies are indeed most concen- 
trated. 

3.3. Diagnostic Diagram 

The presence of an AGN is thought to be best veri- 
fied by the detection of strong high- ionization lin e s such 
as [NeV]14.21/xm or [O iv]25.9^m. iGenzel et al.l (Il998h 
were the first to show that the ratio of high to low exci- 
tation MIR emission lines combined with PAH strength 
could be used to distinguish AGN activity from star- 
formation in ULIRGs . This diagnostic was recently re- 
visited by iDale et all ()2006[ ) for the nuclear and extra- 
nuclear regions of normal star-forming galaxies in the 
SINGS sample observed with the IRS. Dale et al. (2006) 
made use of the [O iv]25.9/im/[Ne ii]12.8/im emission 
line ratios with spline derived EWs of the 6.2/im fea- 
ture (they also proposed an alternative diagnostic us- 
ing the [Si ii]34.8/im/[Ne ii]12.8/im emission line ra- 
tio but [Si ii]34.8^m is beyond the usable range of 
our data). The left panel of Figure [17] shows the 
Dale et al. diagram using the spline derived EWs of 
the 6 .2p,m feature measured by ISargsvan fc WeedmanI 
(|2009l) in the SSGSS sample. The AGN with no de- 



tected [O iv]25.9/xm line is plotted as an upper limit as- 
suming an [O iv]25.9/im/[Ne ii]12.8/im line ratio based 
on the correlation between [Ne ii]12.8/xm/[0 iv]25.9/xm 
and [O iii]A5007/H/3 for other AGNs in Figure [TUl 
We applied a cut in error bars to this plot for clarity 
(Alog([Ne ii]12.8^m/[0 iv]25.9^ni) < 1.5, roughly the 
scale of the y-axis), which excludes 1 AGN, 1 composite 
galaxy and 1 SF galaxy. One other AGN is found with 
no measurable EW. The dotted line represents a variable 
mix of AGN and SF region; the short solid lines perpen- 
dicular to it delineate the AGN region on the left, the 
SF region at the bottom right, and in between a region 
of mixe d classifications w hose physical meaning remains 
unclear ([Dale et al.ll2006[ ). Given the relative homogene- 
ity of our sample (lacking in extreme types), the very 
narrow range of spline EWs for ordinary galaxies and the 
rather large uncertainties in our emission line ratios de- 
rived from low resolution spectra, this diagnostic proves 
of limited use for normal galaxies. Most optically clas- 
sified SF galaxies do fall into the SF corner, but so do 
a few composite galaxies. The rest of the sample shows 
little spread within the mixed region. 

Based on the results of this and the previous sec- 
tions, we revise this diagnostic using the PAHFIT based 
EWs (Eq. [3]) and the correlations between these EWs 
at low wavelength and [N ii]A6583/Hq: (Figure [T5|) on 
the one hand, and between [Ne ii]12.8/im/[0 iv]25.9/im 
and [O iii]A5007/H/3 (Figure on the other 

hand. The right panel of Figure [17] shows the 
[Ne ii]12.8/im/[0 iv]25.9/im emission line ratios against 
the PAHFIT EWs of the 8.6nm feature. Note that 
we inverted the y-axis with respect to the left panel 
(and the traditional Genzel et al. diagram), so that 
the figure becomes a flipped version of the optical BPT 
diagram. The short-dash ed lower line is th e theo- 
retical optical boundary of iKewlev et al.l ()2001fl trans- 
lated into this plane using the correlations between 
EW(8.6Aim) and [N ll]A6583/Ha in Figure [15] and be- 
tween [Ne ii]12.8/im/[0 iv]25.9Aim and [O iii]A5007/H/3 
in Figure [10] for the AGNs and composite galaxies. Its 
analytical form is: 



1.8 



y 



1.51 



-0.88 



(4) 



where x = log(EW(8.6^m)) and y — 
log([Ne ii]12.8Aim/[0 iv]25.9/im). The dotted up- 
per lin e is the empirical boundary of iKaufFmann et ahl 
(|2003D translated using these same correlations for the 
composite and SF galaxies: 



y 



1.10 
x + 0.32 



1.27 



(5) 



As expected from the poorer correlation between 
[Ne ii]12.8Aim/[0 iv]25.9Aim and [O ill] A5007/H/3 for non 
AGNs, this boundary is less meaningful even though it 
does isolate the bulk of the SF galaxies. The long-dashed 
line is an empiral boundary marking the region below 
which we do not find any SF galaxy: 



1.2 



y 



x + 0.8 



0.7 



(6) 



Despite a mixed region of composite and SF galaxies, 
there is a clear sequence from the bottom left to the top 
right of the plots and 3 regions where each optical class is 
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Fig. 17. — Left: The spline derived 6.2/xm PAH equivalent widt hs ISargS^aK fc Weedmaiill2009h against the [O lv]25.9Atm/[Nc n]12.8Atm 
cemission line ratios (a diagnostic diagram originally proposed bv lGenzeletaLI II1998I )). The dotted line represents a variable mix of AGN 
and SF region; the short solid lines perpendic ular to it delineat e the AGN region on the left, the SF region at the bottom right, and 
in betwreen a region of mixed classifications bv lDale et al.l {2006!). We applied a cut in error bars for clarity. This diagram is of limited 
resolving power for normal galaxies. Right: The PAHFIT derived S.Sfim PAH eq uivalent w idths against [Ne ll] 12.8/.tm/[0 iv]25.9/.tm (note 
the reversed y-axis) . This version which resembles a flipped version of the optical IBaldwin~et al.. ( 1981, ) diagram better recovers the optical 
classification. The short-dashed lower line and the dot-dashed upper line are optical boundaries translated into the MIR plane as explained 
in Section 3.3. The long-dashed middle line is an empiral boundary marking the region below which we do not find any optically defined 
SF galaxy. These boundaries are reported in Table |4] The circled galaxy and the lower limit in the SF corner are the two SF galaxies with 
EW lower than the SF limit in Figures [T3l [Til and [TSl (EW(8.6um) <lum). 



uniquely represented. In particular weak AGNs separate 
remarkably well in this diagram. The mixed region may 
in fact be revealing an obscured AGN component in a 
large fraction (> 50%) of the optically defined 'pure' SF 
galaxies. Other dust insensitive AGN diagnostics such as 
X-ray or radio data are necessary to confirm this. Deep 
XMM observations are available only over a small region 
of the Lockman Hole and the FIRST radio limits are 
too bright to reliably test the presence of faint AGNs. 
Indeed we do not expect this hidden AGN contribution 
to be large since none of the SF galaxies falls into the 
AGN corner of the diagram. These objects warrant a 
detailed study beyond the scope of the present paper. 

Alternatively truly 'pure' SF galaxies may be defined 
as lacking the [O iv]25.9/im emission line 25% of our 
SF category). These are not represented except for one, 
which is one of the two SF galaxies with EWs lower than 
the SF limit in Figs [H [H and [l5] (EW(8.6Aim)<l^m). 
The lower limit was calculated by arbitrarily assigning it 
the lowest value of the [OIV] fluxes detected in the sam- 
ple. The other one, which has a detected [OIV] line, is 
circled. These 2 galaxies which would have been mis- 
classified as AGNs based on their EW alone sit well 
into the SF category on this diagram. The AGN with 
no detected [OIV] (plotted as a lower limit) happens to 
have the lowest 8.6/im EW in the sample. A significantly 
larger [Ne ii]12.8/im/[0 iv]25.9/im flux ratio would move 
it into the LINER region of this flipped BPT diagram (al- 
though this particular AGN is not optically classified as a 
LINER). Equations mis] and [6] are reported in Table |4] as 
well as their equivalents for the 6.2/im and 7.7//m PAH 
features. We note that much larger samples, of AGNs 
in particular, are needed to confirm and/or adjust these 



relations. 

4. MIR DUST COMPONENTS AND THE TOTAL 
INFRARED LUMINOSITY 

In this section we investigate how individual dust com- 
ponents emitting in the narrow MIR region trace the 
total dust emission in galaxies, which includes a very 
large FIR component. The definition of the total in- 
frared (TIR) luminosity and the methods used to est i- 
mate it varies in the literature (jTakeuchi et al.l l2005l) . 
In this paper Ltir refers to L(3 — llOO^m) and has 
been derived by fitting the Spitzer photor netric points 
IRAC -URS Blue Peak Up-f-MIPS) with iDraine fc Lil 
2007t) model SEDs" and integrating the best fit SED 
from 3 to 1100/xm. This Ltir is in excellent agreement 
with the 3 — llOOAtm lurn i nosity derived from the pre- 
scription of lDale fc Heloul \2m% for MIPS data (their 
Eq. 4), with a standard deviation of 0.05 dex. This 
shows that the tot al IR luminosity r eally is driven by the 
MIPS points fe.g. iDale eralll2007[ ). We note also that 
integrating the SEDs between 8 and 1000/zm (sometimes 
called the FIR luminosity) would decrease the luminosity 
by ~ 0.04 dex in the present sample. 

Figure [TS] shows the correlations between Ltir and 
LmirI Ltir ratios where Lmir. equals - from top to bot- 
tom - the luminosity of the PAH complexes at 7.7 and 
17/zm, the luminosity of the continuum at 8 and 16/im, 
and the total restframe luminosities in the 8/im IRAC 
band, 16/im IRS band and 24/im MIPS bands. The con- 
tinuum and broadband luminosities are defined as vL^. 
As in all previous figures, SF galaxies are shown as black 

http;/ /www. astro. princeton.edu/~drainc/dust/irem. html 
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Fig. 18. — LmirI Ltir ratios as a function of Ltir where 
Lmir equals - from top to bottom - the luminosity of the PAH 
complexes at 7.7 and 17/im, the luminosity of the continuum at 8 
and 16/jm, and the total restframe luminosities in the 8^im IRAC 
band, 16^m IRS band and 24/im MIPS bands. Symbols are as 
described in Figure [T] The continuum and broadband luminosities 
are defined as vLy. The logarithmic scaling factors k indicated in 
each panel are defined as the median of log(L7ijj{/Ljvi/ii) for the 
SF population alone (green dashed lines). The rms and Pearson 
coefficient r in each panel are also for the SF population alone. 



dots, composite galaxies as pink stars and AGNs as red 
triangles. The logarithmic scaling factors k indicated in 
each panel are defined as the median of \og{LxiR/ Lmir) 
for the SF population alone and is represented by the 
green dashed lines {\og{L m i r / Lt i r) + k = 0). The 
rms and Pearson coefficients of the correlations are also 
quoted for the SF population alone. 

It is striking that galaxies of all types follow the same 
tight, nearly linear correlations between Ltir and the 
broadband luminosities in all 3 Spitzer bands over 2 dex 
in luminosity. This implies that all the galaxies in our 
sample are assigned nearly the same SED from a few 
fim to a thousand /xm and that the FIR component can 
be well predicted from any one broadband luminosity 
in the MIR. This in turn suggests a common heating 
source for the small and large dust grains responsible for 
the M IR and FIR emissions respectively (jRous sel et al.l 
[20nl . The same correlations apply whether this heat- 
ing source is stellar or an AGN. Although this may re- 
sult from the implicit stellar origin of the dust heat- 
ing in the models, the source of ionizing radiation may 
not significantly affect the broad SED, at least for 
weak AGNs. Many attempts have been made to derive 
calibrations b etween Ltir. and sin g le MIR broadband 
luminositi es (ICharv fc Elbazj 120011 lElbaz et al.l 120021: 
Takeuchi et al.l l2005t ISaiina et alT 120061: iBrandl et al.l 
till 



2006t iBavouzet et all 120081: IZhu et al.ll2008ir Our best 
fit slope at 16/xm { Ltir oc ^i6^tfm° '°^) is in good agree- 
ment with that of ICharv fc Elbad ()2001h for the 15^m 
ISO fluxes. At 24/im our correlation for SF galax- 
ies {Ltir oc - ^24um '°^^) is m ore linear than found in 
other studies (iTakeuchi et al.l 12 005: Sai ina et al.l 120061 : 
IZhu et all 120081: IBavouzet et al.l 12008.) but the discrep- 
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Fig. 19. — LmirI Ltir ratios as a function of PAH equivalent 
width at 7.7/im (left panels) and 17/im (right panels) where Lmir 
is defined at the top left of each panel. The symbols and k are 
defined as in Figure [TSl The dotted lines show the expected rela- 
tions when the broadband fiuxes at 8 and IGfira are substituted for 
Ltir in the left and right panels respectively. 



ancy with the firs t thre e calibrations (ITakeuchi et al.l 
l200a ISaiina et al.l I200I iZhu et al.l |2008D disappears 
when including composite galaxies into the fit {Ltir. 
oc ^24^m° °'^)- C)n the other hand our correlation is 
in exc ellent agreement with the iMoustakas fc KennicufrO 
()2006| ) sample (hereafter MK06). 

The PAH and continuum luminosities also correlate re- 
markably tightly and nearly linearly with Ltir, however 
with some distinctions between AGNs and SF galaxies 
and between the hot and cool parts of the spectrum. 
The scatter between Ltir and PAH luminosity for SF 
galaxies is larger for the 17/im PAH feature than for the 
7.7/im PAH feature. AGNs and composite galaxies blend 
with the SF population in the 17/im feature correlation 
whereas they tend to have lower PAH luminosities at 
7.7/ini and stronger 8/.tm continua for the same Ltir- 
The residuals are shown in Figure [T9| as a function of 
the corresponding equivalent widths. The relations be- 
tween these residuals and EWs are of course expected 
since the total flux at 8 and 16/im can be nearly per- 
fectly substituted for Ltir for SF galaxies and AGNs 
alike in the left and right panels respectively (the dotted 
lines show the predicted relations assuming these substi- 
tutions). The most scattered correlation is found with 
the continuum luminosity at 8/im. This may be due to 
larger measurement errors since this continuum is faint 
and/or a stellar contribution unrelated to Ltir- A more 
speculative reason may be that this continuum originates 
from high intensity radiation fields only and is thus un- 
correlated with the cold component of Ltir, unlike the 
PAH emission. 

The scaling factors k are listed in Table |5| for the main 
PAH features and the Spitzer band luminosities. We also 
add to our list of MIR components the peak luminosity of 
the 7.7/im PAH feature, defined as i'L^{7-7fim), as it is a 
more easily measurable quantity at high redshift than the 
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integ rated flux of the PAH feature (jWeedman fc Houc3 
I2OO9I: iSargsvan fc Weedmad I2OO90 . For galaxies with 
EW> 4/xm (most SF galaxies), the median ratio of this 
peak luminosity to the total luminosity of the PAH com- 
plex, h'Li,{7.7^im)/ LpAHi7-7fj,m), is 9.3 ± 0.9 and the 
peak luminosity estimates the total PAH luminosity to 
within ~ 20%. However the overestimate can be as large 
as 50% for other galaxy types in this sample, in particu- 
lar galaxies containing an AGN which may not be easily 
isolated in high redshift samples and may also have much 
smaller EWs leading to yet larger errors. 

For the calibration that shows the strongest deviation 
from linearity in Figure ITSl which is found for the 7.7/im 
PAH luminosity {Ltir oc L^^^f^^'^"^), the linear approx- 
imation log(iT/_R) = log(i_pyiff(7.7^m)) -I- K (where n = 
1.204) recovers Ltir within a factor of 1.5 in this sample. 
For st arbursts and ULIRG starbursts, iRigopoulou et al.l 
dHHI) found a mean \og{LTiR/ Lpah{T IJ-m)) of 2.09 
and 2.26 respectively, conside rably larger than for nor- 
mal galaxies. More recently iLutz et aD ()2003| ) find a 
mean logarythmic ratio of 1.52 for a sample of starburst 
nuclei, closer to our value. Our mean logarythmic ra- 
tio for the 6.2/im feature is 1.5 and 2.0 wi th and with- 
out aperture correction respectively while ISpoon et al.l 
([2004) find a value of 2.4 for a sample of normal and 
starburst nuclei. This ratio is yet higher (3.2) in Galac- 
tic H II regions while highly embedded st ar-forming re- 
gions can lack PAH emission altogether (jPeeters et al.l 
|2004|) . These increased ratios for starburst regions com- 
pared to normal SF galaxies are generally attributed to 
PAH destruction near the site of on-going SF due to in- 
tense radiation fields, making PAHs poor tracers of SF 
(Peeters et al. 2004 and references therein) . The EW de- 
pendence of the \og{LTiR/LpAH{7-7^m)) ratio is clearly 
seen within our sample in the upper left panel of Figure 
[T9I This cautions against the use of a single linear rela- 
tion between PAH luminosity and Ltir for galaxies of 
unknown physical properties. 

However independently of galaxy type we expect to 
find lower values than these studies which all made use 
of interpolation methods to extract the PAH features. 
Using a Lorentzian profile fitting method comparable to 
PAHFIT for a sample of starburst-dominate d LIRGS at 
z ~ 0.5 - 3, iHernan-Caballero et al.l (|2009f ) find mean 
\og{LTiR/LpAH) ratios of 1.92 ± 0.25, 1.42 ± 0.18 and 
1.96 ± 0.27 for the 6.2, 7.7 and 11.3/xm features respec- 
tively. These ratios are 2.6, 1.8, and 1.4 times larger than 
ours respectively, closer than previous studies despite the 
quite different galaxy type considered. The wavelength 
gradient can be explained in the context of selective PAH 
destruction. 

Finally we note that in our sample the total 6.2 — 33/im 
PAH luminosity amounts to ^ 15% of the total IR lu- 
minosity for SF galaxies, ~ 11% for composite galaxies 
and ~ 8% for AGNs. The 7.7/zm feature alone accounts 
for ~ 40% of the total PAH emission. These fractional 
contributions are in good agreement with those found in 
the SINGS sample (S07). 

5. MIR COMPONENTS AND THE STAR-FORMATION 
RATE 

The TIR luminosity is a robust tracer of the SFR for 
very dusty starbursts, whose stellar emission is domi- 
nated by young massive stars and almost entirely ab- 
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Fig. 20. — Lejt: the extinction and r-band aperture corrected 
Ha luminosity against the TIR and 24iim continuum luminosities. 
K is defined as the median of log(Lg'^''/Ljfl) for the SF population 
alone. The rms and slope of the linear regressions (solid lines) are 
also shown for the SF population. The green dashed lines indicate 
equality. The open blue squares and crosses are SINGS data (in- 
tegrated values and galaxies centers respectively); The open green 
squares represent the Moustakas & Kennicutt ( 2O03) (MK06) sam- 
ple. The dotted lines are fits to the MK06 sample. Right: same 
as in the left panels but using the smaller H« aperture corrections 
computed bv iBrinchmann et al.. 12004 ) (see text for details). The 
new correlations (solid lines) are steeper, in better agreement with 
data that do not require aperture corrections. 



sorbed by dust (typically galaxies with depleted PAH 
emission), but for more quiescent and/or less dusty 
galaxies such as those in the present sample, it can in- 
clude a non negligible contribution from dust heated by 
evolved stars ('cirrus emission') as well as miss a non 
negligible fraction of the young stars' emission that i s 
not absorbed by dust (jLonsdale Persson fc Heloulll987D . 
For normal spiral galaxies the contribution of non ioniz- 
ing photons m ay actually domin a te the dust h eating over 
H II regions ()Dwek et al.l I2OOOI: iDwekl I2005D while low 
dust opacity makes these galaxies Ha and UV bright. 
The tight correlations between MIR luminosities and 
Ltir indicate that the same caveats apply from the MIR 
to the FIR (jBoselh et al. 2004). 

5.1. MIR dust and Ha 

Ha emission is a more direct quantifier of young mas- 
sive stars - in the absence of AGN - but inversely it 
must be corrected for the fraction that gets absorbed 
by dust. The SDSS line fluxes are corrected for fore- 
ground (galactic) reddening using IO'Donnel]|(|l99l . The 
correction for intrinsic extinction is usually done us- 
ing the Balmcr decrement and an extinction curve to 
first order, o r more accurately with higher order hy- 
drogen lines (jBrinchmann et 311120041 ). Here we correct 
the SDSS Ha fluxes in the usual simple way using the 
stellar-absorption corrected Ha/H/? ratio and a Galac- 
tic extinction curve. We assumed an intrinsic Ha/H/3 
ratio of 2.86 (case B recombination at electron temper- 
ature Te lOOOOA' and density = 100 cm"^) and 
Ry = A{V)/E{B - V") = 3.1 (the mean value for the 
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diffuse ISM). The Ha attenuations range from 0.4 to 2.3 
mag in the SF galaxy subsample with a median value of 
1.1 mag, meaning that between 10 and 70% of the Ha 
photons do not get reemitted in the IR. 

The SDSS Ha measurements also require fiber aper- 
ture corrections. Here again we apply the usual 
method which consists in scaling the fiber-measured Ha 
fluxes using the r-b and Petrosian-to-fiber flux ratios 
(|Hopkins et al.1 [20031) . The mean value for these ratios 
is 3.5. The left panels of Figure shows the extinction 
and aperture corrected Ha luminosity, L^"' , against the 
TIR and 24^m continuum luminosities (top and bottom 
panel respectively). The rms and slope a of the linear 
regressions (solid lines) are indicated for the SF popula- 
tion alone. The logarithmic scaling factors k indicated in 
each panel are defined as the median of log(Lg'^''/L/fl) 
also for the SF population alone. The green dashed lines 
indicate equality (log(Ly^'') = log(L/fl;) -I- k). Over- 
laid are the MK06 data (open green squares) and SINGS 
data (open blue squares for the integrated measurements, 
cross es for 20"x20" galaxy c enter measurements), taken 
from lKennicutt et all ()2009t ) (hereafter K09). 

Our median Ltir to L^^^ logarithmic ratio of 2.27 ± 
0.2 is in good agreement with the ratio of SFR calibration 
coefRcients derived by Kcnnicutt (1998) for Ha and Ltir^ 
respectively, implying that Ltir (and the MIR compo- 
nents that correlate with it) may be reasonable SFR trac- 
ers in normal SF galaxies after all. This may actually be 
a coincidence due to the fact that the cirrus emission and 
the unattenuated ionizing flux roughly cancel each other 
in massive spiral galaxies (K09 and references therein). 
Our Ltir/L^^^ ratio is also in good agreement with the 
MK06 sample (2.32±0.19). In recent years several groups 
have exploited the capabilities of Spitzer to re- investigate 
the relationship between MIR components and Ha emis- 
sion. Our mean vL^{2AiJLm) to L'^^^ logarithmic ratio 
of 1.31 ± 0.14 is in good agreememt with these (e.; 
iWu et all [20051: iZhu et all [20081: iKennicutt et al.|[200£ 
as is the hig her mean vL ^(2A^im) / L^Jl"^ ratio for compos- 
ite galaxies (IZhu et al.l [20081. However the slope of our 
correlations for SF galaxies tend to be more linear than 
found in these studies (the dotted lines in Figure [501 show 
fits to the MK06 sample). Yet non linearity is expected 
from the positive correlation between attenuation and 
SFR. Given the good agreement between our and the 
MK06 samples in the IR (cf. the Ltir - vL^{2A^m) 
correlation in the previous section), differences in Ha 
measurements must be responsible for the discrepancy 
in slopes. In particular it is possible that aperture cor- 
rections, which are not needed for the MK06 sample, are 
overestimated for all or a fraction of our galaxies. This 
would be the case if SF is more intense at the center 
of the galaxies and/or more a ttenuated, a commo n oc- 
curence in spiral galaxies fe.g. iCalzetti et ai]|2005D . 

As a tes t we consider the smaller aperture corrections 
derived bv lBrinchmann et all (|2004D (hereafter B04) that 
rely on the likelihood distribution of the specific SFR as 
a function of colors. These corrections depend on the 
galaxy colors outside the fiber which are not necessarily 
the same as inside, and are on average 1.6 smaller than 
the r-band corrections for SF galaxies. The right panels 
of Figure [20] shows the same relations as in the left panels 
using these smaller aperture corrections. The new cor- 
relations (solid lines) are indeed steeper while the higher 
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Fig. 21. — The ratio of Ltir to observed Ho luminosity as 
a function of Ha attenuation measured from the Balmer decre- 
ment (SF galaxies only). The top panel assumes conventional 
r-band aperture corrections for Ha, while the bottom panel as- 
sumes the B04 aperture corrections (see text for details). The 
solid lines are best fits to Eq. [7] by K09 for the SINGS-I-MK06 
samples (aTiR = 0.0024 ± 0006). The dashed lines are best fits 
to the SSGSS sample {aTiR = 0.0033 ± 0.0014 in the top panel, 
0.0020±0.0006 in the bottom panel). The smaller aperture correc- 
tions used in the bottom panel significantly improves the fit and 
the agreement between the three samples. 



mean Ltir/ L^^^ and i'L^{24iim) / L^/^^ logarithmic ra- 
tios of 2.47 ± 0.14 and 1.50 dz 0.15 respectively remain 
within the range of the MK06 sample. 

More dramatic is the effect on the relation between 
Ha attenuation and the ratio of Ltir, or other IR lu- 
minosity, to L^^, the 'observed' (aperture-corrected but 
attenuation-uncorrected) Ha luminosity. This relation is 
shown in Figure[2T]for both types of aperture correction. 
K09 modelled the Ha attenuation as: 



Ana = 2.5 log 



1 + aiR 



^IR 
T obs 



(7) 



equivalent to L^^"^ = L^^ + ajp LiR. This en ergy bal- 
ance approach was in troduc ed bv ICalzetti et al . (2007), 
iPrescott et"all ()2007[ ) and IKennicutt et all (|2007[) to 
correct Ha fluxes but has long been used to esti- 
mate UV attenuat ions from the Ltir/ Lpuv ratios (e.g. 
iMeur er et al.lll999l ) . The solid lines in both panels of Fig- 
ure[n]show the best fits by K09 for the SINGS and MK06 
samples {utir — 0.0024 ± 0.006). The dashed lines are 
best fits for the SSGSS sample {aTiR = 0.0033±0.0014 in 
the top panel and 0.0020 ± 0.0006 in the bottom panel). 
The smaller aperture corrections used in the bottom 
panel significantly improves the fit and the agreement 
between the three samples. Unless otherwise stated we 
now assume these corrections in the rest of the paper. 

Figure [22| shows the L^^ -I- uirLir to L^/ZJ" ratios as 
a function of Lff^^ for the TIR, 24//m continuum, 7.7/im 
and 17/im PAH luminosities. The aiR coefficients are 
indicated at the top left of each panel for the SF popula- 
tion. The rms and Pearson coefficients are also indicated 
for the SF population. For the TIR and 24/im luminosi- 
ties, ttTiR = 0.0024 ± 0.0006 and = 0.020 ± 0.005 
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Fig. 22.— L°''^ + ajuLjjt to L^°_^'' ratios as a function of Lg'^'' 
for the TIR, 24/im continuum, 7.7fj,m and 17/im PAH luminosities, 
assuming the B04 aperture corrections for Ha. The am coefficients 
are indicated at the top left of each panel. For the TIR and 24/im 
luminosities, utir = 0.0024 and a24 = 0.020 are best fits to Eq. 
□ for the SINGS+MK07 samples by K09. For the PAH lumosities, 
ay = 0.034 ± 0.012 and an^rn = 0.320 ± 0.159 are best fits to 
Eq.Elfor the SSGSS sample. 



are best fits to Eq. [7] for the SINGS+MK06 samples by 
K09. The combinations of and Ltir or vL^{2Aiim) 
provide a very tight (rms=0.08) and perfectly linear fit 
to the total Ha luminosity for all samples combined over 
5 dex in luminosity, as was shown by K09 for the SINGS 
and MK06 samples. Composite galaxies follow nearly 
the same relation save for 2 over-corrected outliers. Al- 
though more scattered AGNs also follow the SF popula- 
tion. For the PAH luminosities, arj^m = 0.034 ± 0.012 
and ai7^m = 0.320 ± 0.159 are best fits to Eq. [7] for the 
present sample. Here also the combinations of and 
PAH luminosities provide a much improved fit to the 
total Ha luminosity compared to the raw Lpah / L^^^ 
relations (not shown), including for composite galaxies 
and AGNs with the exception of a few outliers, most no- 
tably a composite galaxy with no Ha attenuation and a 
large IR/Ha ratio (#84). 

The same exercise can be performed with similarly 
good results with any other MIR dust components. The 
qtir and a^^nR coefficients for the SSGSS sample are 
listed in Table [5] Note that qmir ^ ^O'^o.tir, using 
the scaling factors k —< \og[LTiRj Lmir) > listed in 
the first column of Table El Although the k factors and 
o-TiR depend on the specific definition of Ltir and on 
the models used to compute it, the umir. coefficients for 
specific dust components or MIR broadband luminosi- 
ties, which are easier to measure than the total IR, are 
independent of these choices. 

As stated above the smaller B04 corrections seem to be 
more appropriate than the usual r-band corrections given 
the agreement with data that do not require aperture cor- 
rections. However they are not trivially calculated (see 
B04 for details of the modeling). More importantly Ha 
is often not easily measurable at all. It is therefore useful 
to provide SFR recipes based on a single IR quantity, or 
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Fig. 23.— Lin/SFRe ratios as a function of SFRe (B04, 
see text for detail) where Lm equals the TIR, 24fim continuum, 
7.7/im and 17/xm PAH luminosities, k is defined as the median of 
log(S F Re / L HI I r) ■ Only SF galaxies for which SFRe is computed 
from the Balmer lines are shown. 



on a combination of IR and UV measurements (see next 
section) as UV is more easily obtained at high redshifts. 
Table[5]lists the mean L^^^ / Lmir. ratios of the SF popu- 
lation for the various MIR components. Keeping in mind 
the non linearities and scatter in the true relations, SFRs 
can be estimated from these approximate Ha luminosi- 
ties using K09's calibration (derived from the latest Star- 
burst99 models and assuming a Kroupa IMF and solar 
metallicity): 

SFRii^ (Mq yr-i) = 5.5 x 10-^2 L'^rr (^^g^ ^-i^^ (g) 

As an example, the SFR derived from the MIPS 
24/im luminosity would be SFR(M0 yr^^) = 6.5 x 
10-^°iyL„{2Afj.m)/LQ consistent with Rieke et al. (2009) 
for galaxies in the range of TIR luminosities of the 
present sample. 

As Eq. [5] was shown by B04 to underestimate the 
SFR of massive galaxies and thus may not be appro- 
priate for this sample or at high 2, we also add to Table 
^SFRe/LMiR calibrations where SFRe is the SFR de- 
rived by these authors as follows: they computed SFR 
likelihood distributions of SF galaxies in the SDSS spec- 
troscopic sample by fitting all stro ng emission lines simul- 
taneously using the ICharlot &: Lo nghetti (2001) models 
and assu ming a Kroupa IMF. Dust was accounted for 
using the ICharlot fc Falll ()2000f ) multicomponent model 
which provides a consistent treatment of the attenuation 
of both continuum and emission line photons. SFR^ 
refers to the medians of these SFR distributions. In this 
model, the Ha attenuation increases with mass while the 
ratio of L^^^ to SFR decreases with mass so that the 
same observed Ha luminosity signals a noticeably higher 
SFR in higher mass galaxies than predicted from Kenni- 
cutt's relation. We refer to B04 for full details. SFRg 
is found to be in good agreement with Eq. [8] for average 
local galaxies but diverges from it for higher mass, higher 
metallicity galaxies such as found in the present sample 
where SFRe is on average twice larger within the SDSS 
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fiber than derived from the Kennicutt relation. How- 
ever the aperture corrections in this study being ~ 1.6 
smaller than derived from the r-band magnitudes for SF 
galaxies, the total SFRe are only ~ 1.3 times larger than 
derived conventionally using the Balmer decrements, r- 
band aperture corrrections and Eq. [H For composite 
galaxies and AGNs, SFRe is not estimated from the 
emission lines which are contaminated by AGN emission, 
but in a statistical way based on the correlation between 
D„(4000) and the specific SFR. We exclude those for 
clarity. 

Figure [23] shows the relations between SFRe and the 
LirI SFRe ratios for the TIR, 24/im continuum, 7.7^m 
and 17/im PAH luminosities. As in previous figures 
the correlation parameters are quoted at the bottom 
right of each panel. These correlations are more scat- 
tered and less linear (higher rms and Pearson coeffi- 
cient) than with Lf^^^. The attenuations underlying 
SFRe being larger than those derived from the Balmer 
decrement for massive galaxies, the SFRe to Ltir ratio: 
SFRe = 3.98 X IQ-'^'^Ltir is very close to that of Ken- 
nicutt et al. (1998) for opaque starburst galaxies (taking 
into account the difference in IMFs). The SFRe/ Lmir 
calibrations are listed in Tabled 

5.2. MIR dust and UV 
Turning now to UV data where dust attenuation is 



an inevitable issue, we once again follow an energy bal- 


ance approach (iMeurer et al.l 


199C; 


Gordon et al. 


200G 


Kone et al.l 120041: iBuat et al.l 


2005 


Cortese et al. 


2008 


Zhu et al.ll2008l: iKennicutt et al.ll2009D. SFRs can 


be es- 



timated from dust corrected FUV luminosities using the 
following calibration by K09 assuming a Kroupa IMF, 
solar metallicity, and adjusted to the GALEX FUV filter 
{Xeff = 1538A). 



SFRfuv (Mq yr- 



. 4.5 X 10-44 (ergs s'^). 

(9) 

where Lff/y = z/L^°''''(1538A) is the dust-corrected 
GALEX FUV luminosity 

Assuming equality with a known SFR estimate (e.g. 
SFRna or SFRe), we derive FUV attenuations as fol- 
lows: 



A 



FUV 



2.5 log 



SFR 



SFRFuv{L%y) 



(10) 



where Lffjy 



vL°J"'{lb?,%k) is the observed FUV lu- 
minosity in ergs s"^. Figure [Ml shows the FUV atten- 
uations of the SF subsample derived from Eq. [TU] as a 
function of Ltir/ L°pfjy (known as the infrared excess 
or IRX) assuming assuming SFR=S'Fi?Ha (Eq. [3 top 
panel) and SFRe (bottom panel). The median FUV at- 
tenuations are 1.9 and 2.8 magnitudes respectively, cor- 
responding to ^ 83 and 92% of the FUV light being 
absorbed by dust (note that assuming conventional r- 
band aperture corrections for Ha yields exactly inter- 
mediate values). The dott ed line is a theoretical re- 
lation by [Buat et all ([2005[): the das hed lines shows a 
model derived bv [Cortese et al.[ ([2008[ ) for galaxies with 
FUV — g = 2.9 corresponding to the mean color of our 
sample (these authors modeled the dependence of the 
IRX/Apijy relation with the age of the underlying stel- 
lar populations, or specific SFR, or color). The solid lines 




FUV/ 

Fig. 24. — The FUV attenuations of the SF population in 
the GALEX FUV band derived from Eq. [10] as a function of 
Ltir/i^L°''"{1530A) (the IRX) assuming SFR = SFRua (Eq. [Sj 
and SFR=S'-F-Ri; (top and bottom panels respectively). The dotted 
line is a theorical relation by Buat et al. (2005); the dashed lines 
shows a model derived by Cortese et al. (2008) for galaxies with 
FUV — g = 2.9 corresponding to t he mean color of our sample; the 
solid lines are best fits to Eq. 1111 
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Fig. 25.— The ratios of FUV to Ha SFRs against the H« SFR: 
Ha is corrected using the Balmer decrement and the B0 4 ap ert ure 
corrections while the FUV is dust corrected using Eqs. llOl and llll 
for the TIR, 24/im continuum, 7.7fim and 17/im PAH luminosities. 



are best fits of the form: 



A 



FUV 



2.5 log 



1 



^IR. 



T obs 
^FUV 



(11) 



equivalent to L'^pjjy 



T obs 
^FUV 



biRLjR, i.e. SFR = 
4.5 X 10-44[L5?^^ + biRLiR], following K09's method. 
Our best fit parameters are bxiR = 0.317 and 0.729 in 
the top and bottom panels respectively. However all 
three models are poor in the bottom panel. FUV at- 
tenuations assuming SFRe are best modeled by a linear 
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Fig. 26. — The ratio of Ne luminosity (defined in Eq. [12} to 
Ltir (a), Lg'^'" (b) and (c) as a function Ltir, L^^'^ and 

SFRe respectively. Only SF galaxies with measured metallicity 
are shown. The lower right panel (d) shows the ratio of the linear 
combination of and Ljve that best fits L^^^'' (see text for 

details) to Lg'^'" against Lg'^''. 



defined as the sum of the 



Fig. 27. — The ratio of H2 luminosity - 
5(0) to 5(2) rotational lines of H2 - to Ltir (a), Lg°^'' (b) and 
SFRf, (c) as a function L^ir, ^ga'^ ^^'^ SFRe respectively. The 
lower right panel (d) shows the ratio of the linear combination of 
Lg'* and Lu., that best fits L??^ 
against Lg°^''. 



(see text for details) to L^"' 



funct ion of log(IRX) or FUV-optical colors (jTrever et alj 
I2007D . Assuming SFRu g (top panel) t he FU V attenua- 
tions are well fit both bv lCortese et al.l ()2008[ ) and by Eq. 



[TTl In this case a linear combination of L°pfjy and Ltib. 
or Lmir recovers SFRna with low scatter as shown in 
Figure [55] using the TIR, continuum, 7.7/xm and 

17/im luminosities. As with Ha in the previous section, 
similarly good corrections can be achieved using other 
MIR components. The hTia. and &A//i? coefficients are 
listed in Table [5] 

5.3. Neon emission lines 



As put forward by iHo fc Ketol (|2007[ ). [Ne ii]12.8^m 
is an excellent tracer of ionizing stars, being an abun- 
dant and dominant species in H 11 regions, quite insen- 
sitive to density, as well as to dust given its long wave- 
length. [Ne lii]15.5/im has similar properties but can be 

the domi nant species in e .g. lo w-mass, low-me tallicity 

galaxies (jCHalloran et alll2006t IWu et al.l 120061 ). Thus 
Ne emission is expected to be directly comparable to the 
dust corrected Ha e mission. Using the CLOUDY code 
()Ferland et al.lll998D . we find that the ionizing flux from 
stars hotter than lOK is best represented by the follow- 
ing weighted linear combination of [Ne ii]12.8/ini and 
[Ne iii]15.5/im also including a metallicity dependence: 



Ha ^ (8.8[Ne II] + 3.5[Ne HI]) x (l/Zf 



(12) 



where Z is the metallicity in solar units. We use the 
right hand side of this equation to define the neon flux 
and luminosity, ^no- 

Fo r the sak e of comparison with the study 

of iHo fc Ketol (|2007( ) who used Ltir as 
SFR estimate, as well as a straight sum of 
[Ne ii]12.8^m and [Ne iii]15.5/.tm, we note that 
our L([Ne ii]12.8^m-|-[Ne iii]15.5^m) to Ltir ratio 
is cons istent with the IRS dataset used by th ese au- 
thors (|0'Halloran et al.l l2006t iWu et al.l 120061) . Our 



I/([Ne ll]12.8^m) to Ltir ratio is larger but this may 
be explained by the large number of low nietallicity 
galaxi es in the samples used, in particular the lWu et al.l 
(|2006f ) dataset which specifically targets low-metallicity 
blue compact dwarf galaxies for which [Ne iii]15.5/im is 
the dominant Ne species (cf. Figure [TT|) . 

The left panels of Figure [26] shows the L^^/Ltir and 
L^c/L^a^ ratios as a function of Ltir (a) and L^^^ (b) 
respectively. Only SF galaxies with measured metallicity 
are represented (85%). Surprisingly Lng behaves much 
like the MIR dust components. It traces fairly linearly 
and tightly the total IR luminosity while we can define 



0.073 ±0.030 using Eq. [Tjsuch that L^^ + aNeL 



Ne 



provides the tightest and most linear correlation with 
^Ha'^i as shown in the lower right panel (d) of Figure [26l 
Likewise we can define b^e — 11.05±5.13 such that 4.5 x 
10~'^'^[L°plfy + bNeLNe] providcs a good fit to SFRna- 
The upper right panel (c) shows the Lnb to SFRe ratio 
against SFR^ , which is significantly more scattered than 
the previous relations as with the MIR dust components. 
This correlation implies the following calibration: 



SFR{Mq yr- 



1.26 X 10-*^ L(Ne) (ergs s"^). (13) 



The and fo^ve coefficient as well as the median 

Ltir/ Lj^c, ^Ha'^/^Ne and SFRe/L^e ratios are reported 
in Table [11 

5.4. Molecular Hydrogen lines 

The rotational H2 lines are fainter than the 
[NeII]12.9^m, [Ne iii]15.5^m and [S iii]18.7^m lines for 
most galaxies in our sample but molecular hydrogen rep- 
resents a significant mass fraction of the ISM in normal 
galaxies. The main excitation source of the rotational 
transitions is expec ted to be FUV radiation from mas- 
sive stars in PDRs (jHoUenbach fc Tielens.1997:. and ref- 
erences therein), therefore these lines also trace SF. The 
first study of warm molecular hydrogen (T ~ 100 — 1000 
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K) in the nu clei of normal, low lu minosity galaxies was 
presented by iRoussel et al.l ()2007D (hereafter ROT) using 
the SINGS sample. A major result of their work is the 
tight correlation between the sum of the S{0) to S{2) ro- 
tational lines (noted F{SO — S2)) and the PAH emission, 
with a ^(5*0 — 52)/PAII ratio insensitive to the intensity 
of the radiation field. This correlation is interpreted as 
supporti ng the origin of H2 exc it ation within PDRs (de- 
fined by Hollenbach fc Tielenj ()1997f) as including the 
neutral ISM illuminated by FUV photons) , with fluores- 
cence as the dominant excitation mechanism. 

Our median logarithmic ratios of L{SO — S2) to the 
TIR, 24/im MIPS band and 7.7/im PAH luminosities for 
the SF population are -3.17 ± 0.19, -2.18 ± 0.23 and 
— 1.95 ± 0.19 respectively. The first two ratios are 1.6 
and 1.8 times larger respectively than those of R07 for 
H II nuclei (taking into account that R07 assumed a fil- 
ter width of 3.1THz instead of vF^ for the 24/im band). 
On the other hand our F{SO - 52)/PAH ratio using the 
stellar component subtracted 8fim IRAC flux instead of 
the PAHFIT extracted feature following R07, is only 1.2 
times higher than that of R07. These differences are 
within uncertainties but the gradients may also reflect 
real differences between H 11 nuclei and disks (warm H2 
more abundant in disks) , as well as support the physical 
link between warm H2 and PAH emissions suggested by 
R07. 

Unlike R07 we do not find significantly higher L{SO — 
S2)/Ltir or L{SO - S'2)/L24 ratios for AGNs (see also 
the top panel of Figure [T0| . This may also be due to the 
much lower AGN contribution when disks are included. 
R07 interpret the higher AGN ratios as an excess of H2 
emission, attributed to additional mechanisms exciting 
H2 molecules in AGNs. The ^(^O - S2)/LpAH{7.7nm) 
ratio does show a significant excess for AGNs, however 
this excess correlates with PAH EWs, suggesting that 
depleted PAHs in AGNs contribute in part to the effect. 
Our dispersion for the ^(5*0 — S2) / Lpah[7 -7 IJ-m) ratio 
is also comparable to the other two while R07 find it 
to be significantly tighter in H 11 nuclei, especially than 
the L{S{) — S2)/ L24, ratio. Our results suggest that local 
complexities are largely washed out on galactic scale and 
that warm molecular hydrogen traces dust in all its forms 
when considering integrated measurements. 

The left panels of Figure [57] show the L{St) - S2) to 
Ltir and L^J^^ ratios as a function of Ltir (a) and L^^^ 
(b). Like Lnc, -^(50 — 5*2) traces reasonably linearly and 
tightly the total IR luminosity while the correlation with 
^Ua' is improved with a linear combination of and 
using Eq. [7] (a//^ — 3.16 ± 2.74) (c). Likewise the 
upper right panel (d) shows the L{SQ—S2) to SFRe ratio 
against SFRe, which implies the following calibration: 

SFR{Mq yr-i) = 6.31 x lO"'*^ L{SO - 5*2) (ergs s^^). 

(14) 

All the coefficients are reported in Table [5] 

6. SUMMARY AND CONCLUSIONS 

We present a MIR spectroscopic survey of 100 'normal' 
galaxies at z ~ 0.1 with the goal of investigating the use 
of mid-infrared PAH features, continuum and emission 
lines as probes of their star-formation and AGN activity. 
Available data include GALEX UV photometry, SDSS 
optical photometry and spectroscopy, and Spitzer near to 
far-infrared photometry. The optical spectroscopic data 



in particular allow us to classify these galaxies into star- 
forming, composite and AGNs, according to the standard 
optical "BPT" diagnotic diagram. The MIR spectra were 
obtained with the low resolution modules of the Spitzer 
IRS and decomposed into unattenua ted features and con - 
tinuum using the PAHFIT code of ISmith et all (|2007t ). 
A notable feature of this decomposition method is to ex- 
tract a much larger PAH contribution (and proportion- 
ally smaller continuum contribution) from the total flux 
compared to standard spline fitting methods which an- 
chor the continuum in the wings of the features where 
non negligible PAH power remain. As a consequence, 
the PAH equivalenth widths are not only larger but ex- 
tend over a considerably larger dynamic range (e.g. the 
equivalent widths of the 6.2 and 7.7/^m PAH features in 
our sample extend to 15 and 32/xm respectively). 

We study the variations of the various MIR spectral 
components as a function of the optically derived age 
(as measured by the 4000A break index), radiation field 
hardness (as measured by the [O iii]A5007/H/3 ratio) and 
metallicity (as measured by [N ii]Af3583/Ha ratio) of the 
galaxies. Systematic trends are found despite the lack 
of extreme objects in the sample, in particular between 
PAH strength at low wavelength and gas phase metal- 
licity, and between the ratio of high to low excitation 
lines (e.g. [O iv]25.9/^m/[Ne ii]12.8/xm) and radiation 
field hardness. These trends confirm earlier results de- 
tected in sources with higher surface brightnesses such as 
ULIRGS, strong AGNs and H 11 nuclei. Our results are 
consistent with the selective destruction in AGN radia- 
tion fields of the smallest PAH molecules efficient at pro- 
ducing the low wavelength PAH features (6.2 to 8.6/im). 
They also suggest that radiation fields harder than those 
in the present sample would also destroy larger PAH 
molecules responsible for the longer wavelength features 
(11.3 to 17/im). Aging galaxies also tend to show weaker 
low wavelength PAH features, consistent with their main 
origin in star-forming regions. 

We revisit the MIR diagnostic diagram of Genzel 
et al. (1998) relating PAH equivalent widths and 
[Ne ii]12.8/im/[0 iv]25.9/im emission line ratios. Based 
on the strongest trends we observed between these mea- 
surements and optical emission line ratios and thanks 
to the extended range of equivalent widths provided by 
PAHFIT, we find this diagnostic to closely resemble the 
optical "BPT" diagram, with a much improved resolving 
power for normal galaxies than previously found based 
on spline derived equivalent widths. A mixed region of 
star-forming and composite galaxies remains, which may 
be revealing obscured AGNs in a large fraction of the 
optically defined 'pure' star-forming galaxies. 

We find tight and nearly linear correlations between 
the total infrared luminosity of star-forming galaxies and 
the luminosity of individual MIR components, includ- 
ing PAH features, continuum, neon emission lines and 
molecular hydrogen lines. This implies that these in- 
dividual MIR components are good gauges of the total 
dust emission on galactic scale despite different spatial 
and physical origins on th e scale of star-forming r egions. 
Following the approach of iKennicutt et al] ()2009l) based 
on energy balance arguments, we show that like the total 
infrared luminosity, these individual components can be 
used to estimate dust attenuation in the UV and in the 
Ha lines. Given the non negligible attenuation in these 
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IR selected galaxies, the correlations between the MIR scaling relations between the various MIR components 
and dust corrected Ha luminosities can also provide first and Ha derived star-formation rates, 
order estimates of the SFR. We thus propose average 
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TABLE 1 

Absorption corrected continuum, PAH and line fluxes derived from the PAHFIT decompositions 
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1.65 ±0, 


.06 


0. 


,49 ±0, 


.02 


5.73 ±0.27 


1, 


.63 ±0.34 


1.24 ±0.25 


6.64 ±0.42 


9 


2, 


.19 ±0, 


.20 


2, 


.26 ±0. 


,06 


6.09 ±0.15 


1.05 ±0, 


.09 


0. 


17 ±0, 


.02 


5.24 ±0.26 


0, 


.16 ±0.20 


0.14 ±0.53 


1.99 ±0.43 


10 


2, 


.11 ±0, 


.07 


1, 


.44 ± 0. 


,09 


3.04 ±0.31 


2.33 ±0, 


.08 


0. 


,38 ±0, 


.02 


5.38 ±0.33 


2, 


.52 ±0.30 


1.02 ±0.46 


3.85 ±0.41 



(a) defined as vFy; 

(b) the sum of H2S(0) to S(2) lines. 



TABLE 2 

Aperture corrections in magnitudes 



ID IRAC Sfira IRSB 16/jm MIPS 2i^im 



1 


0.34 


0.09 


0.06 


2 


0.35 


-0.01 


0.08 


3 


0.29 


0.16 


0.19 


4 


0.86 


-0.04 


-0.08 


5 


0.35 


0.11 


0.16 


6 


0.55 


-0.13 


-0.10 


7 


0.73 


0.02 


0.06 


8 


0.90 


-0.18 


-0.12 


9 


0.54 


0.23 


0.25 


10 


0.46 


0.42 


0.48 



TABLE 3 

PAH EQUIVALENT WIDTHS IN /im 



ID 




6.2/^m 


7.7/im 




8.6/im 




11.3/im 




17^tm 




1 


2, 


.99 ± 0. 


,38 


8.50 ±0.73 


1, 


.80 ±0. 


,22 


2, 


.86 ±0.28 


1, 


.19 ±0, 


.12 


2 


4, 


.82 ± 1. 


39 


11.30 ±2.06 


2, 


.44 ± 0. 


,55 


3, 


.07 ±0.47 


3, 


.00 ±0, 


.39 


3 


2, 


.91 ±0. 


,27 


9.15 ±0.77 


1, 


.90 ± 0. 


,18 


4, 


.79 ±0.31 


1, 


.71 ±0, 


.29 


4 


2, 


.13 ±0. 


,92 


6.37 ± 1.60 


1, 


.74 ±0. 


,59 


2, 


.08 ±0.53 


1, 


.81 ±0, 


.54 


5 


2, 


.59 ±0. 


,57 


12.92 ± 1.10 


3, 


.00 ± 0. 


,31 


4, 


.53 ±0.33 


2, 


.70 ±0, 


.43 


6 


2, 


.00 ± 0. 


,34 


5.67 ±0.68 


1, 


.41 ±0. 


,29 


1, 


.99 ±0.36 


1, 


.91 ±0, 


.46 


7 


2, 


.55 ±0. 


,44 


10.95 ±0.89 


2, 


.13 ±0. 


,16 


2, 


.37 ±0.13 


1, 


.64 ± 0, 


.17 


8 


1, 


.74 ± 0. 


17 


6.48 ± 0.30 


1, 


.73 ± 0. 


,07 


2, 


.09 ± 0.06 


1, 


.44 ± 0, 


.07 


9 


0, 


.93 ±0. 


27 


3.70 ± 0.49 


0, 


.67 ±0. 


10 


2, 


.52 ±0.20 


1, 


.09 ± 0, 


.13 


10 


2, 


.17 ±0. 


,15 


8.48 ± 0.45 


1, 


.60 ±0. 


,11 


4. 


.27 ±0.20 


4, 


.17 ±0, 


.46 



TABLE 4 

Analytical boundaries in the plane of PAH equivalent widths [/^m] versus [Ne ii]12.8/iM/[0 iv]25.9/iM("'. 





6.2/im 


7.7/tm 


8.6^m 




Eq. [4] Eq. [5] Eq. |6] 


Eq. H Eq. \5\ Eq. \E\ 


Eq. H Eq. [5] Eq. |6] 


Cl 


2.36 1.96 1.00 


2.21 1.36 1.30 


1.84 1.10 1.20 


C2 


1.84 0.64 0.90 


1.16 -0.21 0.55 


1.51 0.32 0.80 




-0.87 -1.30 -0.35 


-0.88 -1.30 -0.45 


-0.88 -1.27 -0.70 



(a) log([NeII]/[OIV]) = [ci/(log(EW) + ca)] + eg. 
Eqs |4l [5] and |6] are defined in Section 3.3. 
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TABLE 5 

Median ratios for the star-forming population 



\og[LTIR/LMIR] ^Og\L']^7/LMIRT^ log[L"^°;7LM/fl]"'^ log[5Ffle/LM/fl]'"' g^f^j b\^. 



TIR 







± 







-2.266 


±0.159 


-2.471 ±0. 


144 


-43, 


.461 


± 


0. 


183 


0, 


.002 


± 


0. 


,001 


0.318 ± 0.124 


PAH ^.^^lrI^ 


1, 


,204 


± 


0, 


,087 


-1.091 


±0.181 


-1.242 ±0. 


154 


-42, 


.248 


± 


0. 


,213 


0, 


.034 


± 


0. 


,012 


5.154 ±2.299 


peak 7.7fj,m 


0, 


,242 


± 


0, 


,086 


-2.054 


±0.175 


-2.231 ±0. 


140 


-43, 


.221 


± 


0. 


,205 


0, 


.004 


± 


0. 


,001 


0.539 ± 0.229 


PAH 11.3/im 


1, 


,851 


± 


0, 


.071 


-0.458 


±0.174 


-0.611 ±0. 


146 


-41, 


.605 


± 


0. 


183 


0, 


.143 


± 


0. 


,053 


21.940 ±9.517 


PAH 17/im 


2, 


162 


± 


0, 


.121 


-0.101 


± 0.203 


-0.292 ± 0. 


192 


-41, 


.263 


± 


0. 


,212 


0, 


.320 


± 


0. 


160 


46.986 ± 29.721 


IRAC 8/im 


0, 


,594 


± 


0, 


.075 


-1.719 


±0.177 


-1.896 ±0. 


142 


-42, 


.871 


± 


0. 


195 


0, 


.008 


± 


0. 


,003 


1.214 ±0.514 


IRSB 16/im 


1, 


,040 


± 


0, 


.058 


-1.253 


±0.157 


-1.445 ±0. 


140 


-42, 


.409 


± 


0. 


185 


0, 


.022 


± 


0. 


,007 


3.544 ± 1.381 


MIPS 24^lm 


0, 


,959 


± 


0, 


.096 


-1.311 


±0.142 


-1.505 ±0. 


151 


-42, 


.508 


± 


0. 


191 


0, 


.019 


± 


0. 


,007 


2.763 ± 1.253 




1, 


,537 


± 


0, 


.183 


-0.738 


±0.216 


-0.930 ±0. 


199 


-41, 


.901 


± 


0. 


,298 


0, 


.073 


± 


0. 


,031 


11.258 ±5.228 


H2S{0)-S{2)(9) 


3, 


172 


± 


0, 


.191 


0.888 ±0.255 


0.706 ±0.239 


-40, 


.272 


± 


0. 


,267 


3, 


.161 


± 


2. 


,641 


411.216 ±486.911 



(a) using r-band aperture corrections 

(b) using B04 aperture corrections (see text for detail). Note that the relations in this case are markedly non-linear (right panel of Figure [20t . 

(c) SFRe [M0/yr] from B04 (see text for detail); Lmih [erg/s] 

(d) Eq. m 

(e) Eq. \m 

if) Eq. \m 

(g) the sum of H2S(0) to S(2) lines 



